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Chapter  1 

INTRODUCTION 

STATEMENT  OF  THE  PROBLEM 


GENERAL 

Historically,  geotechnical  failures  caused  by  earthquakes  have  been  a  serious  problem 
which  is  not  fully  understood.  Stories  of  liquefied  soil  bubbling  out  of  the  ground,  both  during 
and  after  ground  motion  has  cea.sed,  have  been  noted  in  numerous  accounts  of  hi.storic 
earthquakes.  Because  of  the  extensive  loss  of  life  and  tremendous  amounts  of  damage  that  these 
earthquakes  and  related  ground  failures  have  caused  during  this  century,  mans  investigators 
have  now  started  studying  these  natural  phenomena.  In  the  last  10-15  years  the  level  of 
understanding  the  mechanisms  of  liquefaction  have  moved  from  an  initial  (pialitative  stage,  to 
recently  a  much  more  detailed  qualitative  and  now  a  (juantitative  level.  The  ultimate  goal  of 
every  r  jrcher  in  the  field  is  the  hope  of  being  able  to  predict,  and  even  prevent  conditions  of 
liquefaction  both  accurately  and  cost  effectively.  There  are  still,  unfortunately,  many  ()uestions 
pertaining  to  our  present  level  of  knowledge  that  must  be  addre.ssed  before  we  may  eliminate 
the  damage  potential  of  this  phenomenon. 

The  basic  concepts  of  liquefaction  were  first  illustrated  by  Ca.segrande  (1936)  in  his 
studies  on  slope  stability.  Casegrande's  work  dealt  mainly  with  volume  change,  and  the 
corresponding  pore  pre.s.sure  changes  during  the  application  of  a  .static  shearing  stress.  In  recent 
years  cyclic  stress,  like  those  created  by  earthquakes  have  been  found  to  cause  similar  volume 
changes  with  resulting  pore  pressure  phenomenon.  This  study  deals  with  the  potential  for 
increa.sed  pore  pressure,  and  the  subsequent  threat  of  liquefaction,  created  by  cyclic  stres.ses 
from  earthquake  generated  ground  motion  in  a  soil  profile. 
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DEFINITION  OF  LIQUEFACTION 

The  term  liquefaction  has  been  used  and  defined  by  different  investigators  over  the  years. 
There  is  also  much  controversy  over  how  the  term  should  be  used  to  describe  different 
phenomena.  The  term  liquefaction,  as  it  is  used  in  this  .study,  will  refer  to  the  changing  of  a 
.soil  from  a  solid  state  to  a  liquefied  state  due  to  the  build-up  of  excess  hydrostatic  pore 
pressure.  This  change  of  state  is  further  defined  as  a  result  of  a  loss  of  shear  strength  due  to 
cyclic  loading  involving  the  reversal  of  shear  stresses  during  dynamic  loading,  such  as  an 
earthquake.  VV'hen  a  saturated  fine-to-medium  grained  cohesionless  soil  mass  is  subjected  to 
ground  vibration,  it  tends  to  compact  and  decrease  in  volume.  If  drainage  of  the  pore  spaces  is 
unable  to  occur,  the  tendency  to  decrease  in  volume  is  met  with  hydrostatic  resistance  and  the 
result  is  an  increase  in  pore-water  pressure.  When  this  excess  pore-water  pressure  builds  to  the 
point  at  which  it  is  equal  to  the  overburden  pressure,  the  effective  stress  becomes  zero,  the 
material  looses  its  shear  strength,  and  the  resulting  state  is  one  of  a  dense  fluid.  At  this  point 
the  material  is  said  to  have  experienced  initial  liquefaction.  Seed  (1976)  describes  the  build-up 
of  excess  pore  pressure  as  the  trade  off  between  a  tendency  for  a  volume  decrease,  due  to  cyclic 
loading,  and  a  rebound  of  the  soil  structure  due  to  the  load  being  transferred  from  the  soil 
.structure  to  the  pore  water.  The  difference  between  soil  structure  rebound  and  the  tendency 
for  the  volume  of  the  soil  mass  to  compact  results  in  the  excess  pore  pressure. 

This  description  of  liquefaction  appears  in  equation  form  as  follows: 

7ff  =  (j'ff  tan  <t>  =  (cTff  -  pi.)  tan  <b 

where 

<6  =  angle  of  internal  friction 

rff  =  shear  stress  on  the  failure  plane  at  failure 

aff  =  normal  stress  on  the  failure  plane  at  failure 

H  =  pore-water  pressure 

H  refers  to  n  static  +  n  excess  build-up. 

Then  as  approaches  off  due  to  dynamic  loading. 

Tff  approaches  0.  which  responds  as  a  fluid. 
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CASE  HISTORIES 

The  sudden  interest  in  expanding  our  knowledge  of  the  liquefaction  phenomenon  can  be 
credited  to  a  select  few  major  earthquakes  occurring  in  the  last  15-20  years.  The  most  dramatic 
of  these  events  occurred  on  June  16,  1964.  On  this  da\'  Niigata.  Japan,  was  hit  by  a  shock 
which  registered  7.3  on  the  Richter  scale  with  its  epicenter  located  35  miles  north  of  the  city 
(Seed  and  Idri.ss,  1967).  During  the  event,  the  extcn.sive  liquefaction  of  the  sand  deposits 
resulted  in  major  damage  to  buildings,  bridges,  highways,  utilities,  dock  areas,  oil  refineries, 
and  railroads.  The  most  dramatic  case  of  settlement  and  tilting  was  undoubted!)  that 
experienced  during  a  major  bearing  capacity  failure  by  apartment  buildings  that  tilted  up  to 
80°  from  vertical.  The  occupants  of  the  complexes  were  then  able  to  walk  down  the  face  of  the 
structure  to  safety  (Seed,  1969).  During  the  same  event  a  sewage  treatment  tank  floated  to  the 
surface  because  the  buoyant  forces  lifting  the  tank  were  not  resisted  by  the  soil  uhich  had 
liquefied  In  another  part  of  the  city  a  flowing  spring  3-feet  high  developed  in  an  urbanized 
area  as  a  result  of  liquefaction  and  a  release  of  the  excess  pore-water  pressures  to  the  surface. 

Another  illustrative  example  of  extensive  damage  created  by  liquefaction  is  the  1964 
Alaskan  earthquake.  This  large  magnitude  earthquake  (8.4)  struck  .South  Central  Alaska, 
disrupting  many  highways  and  utilities  and  severely  damaging  structures.  Voud  (1978) 
estimates  that  ground  failure  caused  60%  of  the  estimated  $300  million  (1964  value)  total 
damage.  Five  major  landslides  caused  about  $50  million  of  the  $85  million  damage  to 
nonmilitary  facilities  in  the  city  of  Anchorage.  Lateral  spread  failures  damaged  highwas-  and 
railway  grades  and  bridges,  requiring  about  $50  million  in  repairs,  and  severely  disrupting 
these  lifelines.  Flow  failures  in  three  Alaskan  coastal  communities  carried  away  ports,  adjacent 
warehouses,  and  transportation  facilities  that  cost  approximately  $15  million.  In  each  of  these 
ground  failure  modes  the  probable  cause  can  be  traced  to  liquefaction.  Ross,  Seed  and 
Migliaccio  (1969)  also  reviewed  a  number  of  highway  bridge  facilities  and  indicated  that 
liquefaction  of  the  support  soils  during  and  shortly  after  the  earthquake  cau.sed  the  failures. 
Turnagain  Heights,  an  upper-middle  class  residential  area  in  Anchorage,  slid  into  the  sea  as  a 
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result  of  the  liquefaction  of  thin  (2-3  inches)  sand  and  silt  lenses  within  the  sensitive  Bootleg(jer 
Cove  clay. 

In  1971  an  earthquake  (M  =  6.6)  struck  the  San  Fernando  Valley  of  California  causing 
millions  of  dollars  of  damage  (Youd.  1971).  A  large  portion  of  this  damage  was  attributed  to 
liquefaction.  The  near  catastrophic  failure  of  the  lower  San  Fernando  Dam  has  been  attributed 
to  liquefaction  of  the  embankment  soil  (Seed,  et  al,  1975).  The  potential  for  a  cata.strophic  loss 
of  life  was  ever-present  due  to  the  dam’s  location  above  a  large  population  center.  Lew, 
Levendecker  and  Dikkers  (1971)  describe  the  extensive  damage  suffered  by  a  juvenile  hall  due 
to  ground  cracking  and  ground  spreading  as  a  result  of  liquefaction  in  the  foundation  soils.  The 
Jo.seph  Jen.sen  Filtration  Plant  also  suffered  structural  damage  due  to  slope  failure,  ground 
surface  cracking,  and  differential  building  settlement  (Dixon  and  Burke,  1973). 

Youd  and  Hoose  (1976)  de.scribe  various  ground  failures  that  occurred  during  the  1906  San 
Francisco  earthquake  that  have  now  been  attributed  to  liquefaction  of  the  soils.  Estimates  of 
85%  of  the  damage  to  San  Francisco  was  caused  by  the  fire  that  followed  the  earthtjuake. 
Much  of  the  fire  damage,  however,  could  have  been  prevented  if  several  main  water  lines  into 
the  city  had  not  been  severed  bv  a  liciuefaction-induced  ground  failure. 

TECHNOLOGICAL  ADVANCES 

As  pointed  out  by  the  previous  examples,  liquefaction-induced  ground  failures  have  been 
directly  or  indirectly  responsible  for  many  expensive  and/or  fatal  problems  a.ssociated  with 
earthquakes.  As  a  result,  renewed  intere.st  among  noted  investigators  has  been  initiated  in  the 
field  of  liquefaction  potential. 

Laboratory  testing  procedures,  such  as  the  cyclic  triaxial  test,  cyclic  simple  shear  test,  and 
the  shaking  table  test  have  greatly  enhanced  the  formulation  of  basic  theories.  Data  from  these 
tests,  although  sometimes  ambiguous  as  to  their  interpretation,  can  enable  investigators  to 
determine  the  causes  of  liquefaction,  the  factors  influencing  the  mechanism,  and  procedures  of 
analysis  and  design  for  the  control  of  the  liquefaction  problem. 


As  in  any  science,  initial  theories  are  proposed  and  tested  against  case  histories  and  field 
observations.  Further  problems  are  identified  and  a  back-anahsis  approach  ma\  be 
implemented  to  formulate  empirical  methods  of  analysis.  A  new  source  of  data,  which  may 
prove  very  important  in  the  future  is  the  available  information  from  the  Chinese  records  on 
earthquakes.  It  is  estimated  that  these  records  may  extend  further  back  in  histors  than  an\ 
other  known  earthquake  records.  The  combination  of  these  records  with  our  present  da\ 
earthquake  records,  and  subsequent  liquefaction  observations  creates  a  substantial  collection  of 
information  which  will  be  available  for  analysis  and  comparison. 

Although  much  work  has  been  done  in  the  prediction  of  liquefaction  potentials,  it  must  be 
noted  that  a  universally  accepted  method  of  analysis  has  not  yet  been  developed.  Any  and  all 
information  that  can  be  collected  and  used  to  find  an  acceptable  approach  to  the  problem 
would  be  a  major  contribution  toward  maintaining  the  integrity  of  structural  foundations 
during  earthquakes. 


PURPOSE  OF  THIS  STUDY 

The  purpose  of  this  study  is  to  develop  a  useful  guide  for  evaluating  the  possibility  of 
liquefaction  and  the  subsequent  potential  for  bearing  capacity  failures  and/or  differential 
settlement  of  dams,  buildings,  and  other  structures  proposed  by  the  United  States  Navy  at  their 
Weapons  Research  Center  in  China  Lake.  California.  This  study  is  meant  to  be  a  contribution 
to  the  work  already  being  undertaken  by  the  Geological  Sciences  group  in  an  effort  to  respond 
to  “Project  2000;"  a  generalized  site  investigation  of  the  base  area  outlining  the  Navy's 
potential  for  expansion  through  the  year  2000.  This  potential  map  will  be  useful  as  a 
preliminary  tool  in  identifying  general  areas  of  potential  liquefaction  and  suggested  design  life 
criteria  for  the  proposed  structures  in  regards  to  dynamic  stability  of  the  soil.  At  the  time  of  a 
site’s  proposal,  depending  on  the  project,  a  more  detailed  site  investigation  should  be 
performed  for  any  site  where  liquefaction  is  considered  possible.  Threshold  accelerations,  (ap){, 
and  a  subsequent  potential  ratio  (3p/(  ap)j,  design  acceleration  divided  by  threshold 


acceleration)  are  (luantitative  measures  for  comparison,  and  should  not  be  considered  finalized 
design  values  for  major  structures. 


The  Naval  Weapons  Center  is  located  in  a  very  seismically  active  area.  The  area  has  been 
classified  as  a  Zone  4  region,  the  highest  cla.s.sification,  by  the  1976  Uniform  Building  Code's 
.seismic  zone  cla.ssification  map  (Figure  1.1).  This  zone  is  outlined  by  the  UBC  as  an  area  of 
major  damage,  capable  of  an  Intensity  VIII  or  greater  on  the  modified  .Mercalic  scale.  In 
addition,  this  area  is  said  to  be  capable  of  further  danger  potential  due  to  its  proximits  to 
major  fault  sy.stems.  Roquemore  (1981),  staff  research  geologist  at  the  Naval  Weapons  Center. 
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VII  Of  THE  M  M  *  SCALE 
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AND  HIGHER  OF  THE  M  M  *  SCALE 
ZONE  4  .  THOSE  AREAS  WITHIN  ZONE  NO  3  DETERMINED  BV 

THE  PROXIMITY  TO  CERTAIN  MAJOR  FAULT  SYSTEMS 

•MODIFIED  MERCALIC  INTENSITY  SCALE  OF  1931 


FIGURE  I.l.  1976  Uniform  Building  Code’s  Seismic  Zone  Cla.ssification  Map. 


is  currently  engaged  in  a  detailed  fault  mapping  program  which  will  locate  all  major  fault 


.sy.stems  in  the  base  vicinity  and  analyze  their  potential  for  activity.  His  on-going  work  will  be  a 


constant  source  of  updated  input  for  the  computer  generated  probabilistic  peak  accelerations  to 
be  used  in  the  future  prediction  of  liquefaction  potential. 

Numerous  base  facilities  constructed  in  the  last  10-15  years  have  been  subjected  to 
differential  settlement  resulting  in  cracked  walls  and  undulating  floors.  It  is  realized  that 
fluctuations  in  the  already  shallow  water  table  may  alter  effective  stresses  causing  consolidation 
and/or  settlement.  However,  .selected  fault  trenches  analyzed  by  both  Roquemore  (1981)  and 
St.  Amand  (1981)  have  illustrated  sand  blows  or  mud  volcanoes  at  depth,  indicating 
liciuefaction  and  potential  foundation  instability.  (See  Figure  1.2.) 


FIGURE  1.2.  Mud  Volcano  in  Little  Lake  Fault  Trench  Adjacent  to  Study  Area. 

Because  of  the  high  seismicity  of  the  area,  a  high  groundwater  table,  and  the  fact  that 
many  soils  at  China  Lake  are  highly  susceptible  to  liquefaction,  it  is  very  important  that 
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potential  locations  for  liquefaction  be  earmarked.  Sudden  and  possibly  catastrophic  failures  of 
structures,  toxic  waste  settling  ponds,  and  soil  embankments  could  possibh'  be  presented  if 
licpiefaction  susceptible  areas  are  outlined.  Given  an  understanding  that  potential  sites  could 
liquefy,  more  specific  design  parameters  could  be  implemented  to  correct  problems  that  exist  at 
a  site,  or  a  problem  site  may  be  avoided  completely. 

The  creation  of  “Project  2000“  indicates  the  Navy’s  anticipated  expansion  of  the  facilits'  at 
China  Lake.  It  is  the  goal  of  this  study  to  aid  the  planners  and  deselopers  in  both  their  site 
.selections  and  the  level  of  insestigation  and  design  that  a  specific  site  might  require. 

STUDY  AREA 

The  Naval  VV'eapons  Center  (NWC)  at  China  Lake  is  a  federal  research  facility  employing 
approximately  3000  civilians  under  the  Department  of  Defense.  Significant  projects  in  the  past 
have  included  the  development  of  the  Sidewinder  missile,  numerous  air-to-air,  air-to-ground, 
and  ground-to-ground  guidance  .systems,  and  a  very  large  collection  of  other  weaponry 
research  projects  either  too  numerous  to  list  or  the  publication  of  them  would  be  classified  for 
security  rea.sons. 

A  survey  of  the  structures  currently  on  site  encompases  3  very  large  multistoried 
laboratories,  numerous  warehouse  facilities  and  operations  and  control  buildings.  Support 
facilities,  including  housing  and  special  services  were  constructed  for  a  capacity  population  of 
nearly  10,000.  Armitage  Air  Field  is  the  base  for  a  moderately  sized  attack  squadron  as  well  as 
the  host  for  numerous  other  military  aircraft  temporarily  assigned  to  maneuvers  and  tests  run 
in  conjunction  with  NWC.  Also  located  on  .site  is  the  Supersonic  Naval  Ordnance  Research 
Track  (SNORT)  which  is  a  5-mile  track  utilized  for  high  speed  rocket  and  sled  launches, 
ballistic  tests  and  many  other  sensitive  operations  rec}uiring  rigid  te.st  conditions  and  a 
dependable  soil  foundation  for  the  track. 

The  base  and  range  area  is  located  approximately  100  miles  due  ea.st  of  Bakersfield, 
California.  The  facility  stretches  some  1800  square  miles  across  the  high  desert  region  at  the 
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so\ithmi  fncl  of  the  Sierra  Nevada  mountain  range.  Due  to  onl\  a  locally  high  groundwater 
table.  howe\er,  this  stud\-  will  be  restricted  to  a  260  square  mile  region  located  at  the  southern 
c'nd  of  the  main  facilits'  area. 


METHODOLOGY 

The  methodology  utilized  in  this  study  is  outlined  b\-  the  flowchart  illustrated  in 
Figure  1.3.  This  method  will  combine  concepts  introduced  by  Youd  and  Perkins  (1978)  with 
tho.se  of  Dobry  and  others  (1979).  In  combining  the  two,  an  attempt  will  be  made  to  utilize 
Youd  and  Perkins’  ability  to  qualitatively  and  empirically  predict  ground  failures  due  to 
liquefaction  over  a  large  area,  and  Dobry’s  quantitative  approach  to  the  prediction  of 
liquefaction  potential  by  the  generation  of  excess  pore  pressure  in  the  soil  profile. 

First,  a  detailed  geologic  map  of  the  study  area  is  used.  Particular  consideration  is  given  to 
the  composition  of  each  unit,  such  as  the  soil  type,  method  of  deposition,  age  of  the  deposition, 
and  any  particular  cementation  noted.  Other  geotechnical  factors  affecting  soil  susceptibilits  to 
liquefaction  are  not  considered  at  this  stage.  The  result  of  this  input  will  enable  the  deposit  as 
a  unit  to  be  cla.ssified  as  to  low,  moderate,  or  high  susceptibility. 

Depths  to  the  groundwater  table  were  then  collected  from  existing  well  sites  on  ba.se.  In 
areas  lacking  such  information,  a  literature  search  was  carried  out  to  determine  previous  levels 
and  current  rates  of  table  fluctuations.  These  depths  and  well  locations  were  digitized  for 
determination  of  preci.se  latitudinal  and  longitudinal  coordinates  before  being  formatted  for 
computer  input.  Two  programs  were  utilized  by  the  computer  in  the  generation  of  a  depth  to 
groundwater  contour  map.  The  first  program  assembled  all  input  data  points  and  created  a 
rectangular  grid  .system  with  calculated  depths  at  all  intersections  or  grid  corners.  The  original 
data  points  were  then  disregarded.  The  .second  program  was  a  graphic  tool  that  connected  all 
grid  points  of  equal  depths  with  contour  lines.  The  contour  lines  ranged  from  water  levels 
located  above  ground  surface  to  depths  of  40  feet,  at  2-foot  intervals.  The  accuracy  of  these 
griding  and  contouring  programs  were  then  verified  by  the  proximity  of  each  original  input 


FIGURE  1.3.  Methodology  Utilized  in  the  Study  of 
Liquefaction. 


II 


data  point  with  the  nearest  contour  line.  This  graphic  procedure  was  not  accepted  until  each 
contour  line  accurately  repre,sented  the  originally  known  input  data  points. 

All  known  boring  logs  were  then  collected  and  the  down-hole  profile  was  analyzed  for 
potentially  liquefiable  soils.  All  borings  revealing  a  lack  of  saturated  sand  or  silt  within  the 
upper  40  feet  enabled  the  locality  to  be  classified  as  non-liquefiable. 

The  depth  to  groundwater  contour  map  was  then  superimposed  over  the  geologic  map 
with  the  units  cla.ssified  as  "high”  and  “moderate"  highlighted.  These  2  cla.ssifications,  when 
falling  in  the  range  of  depths  to  water  less  than  or  equal  to  40  feet  were  then  recorded  as 
potentially  prime  liquefaction  areas.  It  is  these  areas  that  were  further  analyzed  in  a  more 
quantitative  manner. 

Dobry  and  others  (1979)  have  proposed  a  profile  stiffness  method  for  liquefaction 
potential.  This  method  utilizes  the  empirical  and  experimental  findings  of  many  authors  to 
develop  the  threshold  peak  ground  surface  accelerations.  (ap)(,  at  which  the  soil’s  cyclic  strain 
is  equal  to  a  strain  level  which  has  been  found  by  many  authors  to  be  the  level  abo\e 

which  excess  pore  pressures  develop.  This  will  be  discussed  in  detail  in  Chapter  4.  Field  data 
required  for  calculation  of  the  threshold  peak  ground  surface  acceleration  includes  the  shear 
wave  velocity  and  the  mass  density  of  the  soil.  An  ES-1200  twelve  channel  signal  enhancement 
sei.smograph  was  utilized  in  the  determination  of  the  .shear  wave  velocities  from  surface 
refraction.  Available  borings  and/or  a  sand  cone  apparatus  were  utilized  in  the  determination 
of  the  .soil's  mass  density.  Both  of  these  values  were  averaged  over  the  entire  geologic  unit, 
thereby  defining  the  threshold  acceleration  for  a  single  geologic  unit  as  a  function  of  only  the 
depth  to  groundwater  and  the  resultant  total  overburden  stress.  These  threshold  acceleration 
values  were  calculated  at  .sites  on  half  mile  centers  in  a  grid  system  set  up  over  the  entire  studs- 
area. 

A  .second  grid  .system,  identical  to  the  first,  was  incorporated  to  analyze  the  seismicity  of 
the  study  area.  A  computer  program  called  EQRISK  (McGuire,  1976)  was  then  implemented  to 
predict  peak  ground  surface  accelerations  on  this  grid.  EQRISK  generates  a  probabilistic  peak 
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ground  surface  acceleration  knowing  the  fault  characteristics  in  the  area,  their  location,  and 
the  region’s  historic  background  seismicity. 

Values  of  the  calculated  threshold  peak  ground  surface  accelerations.  (ap)(.  and  the  design 
peak  ground  surface  accelerations,  ap,  were  compared  for  a  100  year  return  interval.  The 
result  was  a  regional  map  illustrating  the  potential  and/or  the  relative  ease  of  pore  pressure 
increases  and  the  po.ssible  liquefaction  at  any  specific  area. 
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Chapter  2 

REVIEW  OF  THE  LITERATURE 

GENERAL 

The  realization  that  liquefaction  is  seismically  induced  is  relati\ely  new.  In  1930 
Casegrande  introduced  the  term  “liquefaction"  as  it  related  to  slowly  applied  loads.  The 
present  concern  over  cyclic  loading  causing  liquefaction  only  began  in  the  early  1960’.s  as  a 
direct  result  of  the  case  histories  previously  described.  Since  that  time  numerous  investigators 
have  supplied  answers  to  many  questions  as  to  what  liquefaction  is,  what  causes  it.  what 
factors  influence  its  development,  and  how  to  alleviate  or  design  for  the  problem.  In  the  past 
10-15  t'ears  there  has  been  much  literature  generated  regarding  laboratory  studies,  case 
histories,  and  methods  of  analysis  of  the  liquefaction  problem.  The  purpose  of  this  chapter  is 
not  to  review  all  the  details  of  these  studies,  or  even  to  attempt  a  complete  list  of  what  has 
been  accompli.shed.  This  chapter  will  describe  liquefaction  and  present  various  methods  of 
analysis  that  are  now  being  utilized,  as  well  as  introduce  findings  that  have  recently  been 
developed. 

DEFINITIONS  OF  LIQUEFACTION 

Because  of  the  numerous  independent  studies  that  have  been  performed  on  liquefaction, 
there  have  been  some  slight  differences  in  definitions  of  terms.  Youd  (1975)  observes  two  basic 
components  to  most  definitions  of  liquefaction;  a  strength  loss  criteria  and  a  flow  deformation 
criteria.  He  indicates  the  attempt  to  combine  these  two  distinct  phenomena  into  a  single 
universal  definition  has  resulted  in  controversy.  Some  of  the  definitions  that  Youd  uses  to 
illustrate  his  point  are: 
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"The  sudden  decrease  of  shearing  resistance  of  a  quicksand  from  its  normal 
value  to  almost  zero  without  the  aid  of  seepage  pressure."  (Terzaghi  and 
Peck,  1948) 

"The  sudden  large  decrease  of  shearing  resistance  of  a  cohesionless  soil, 
cau.sed  by  a  collapse  of  the  structure  by  shock  or  strain,  and  associated  with 
a  sudden  but  temporary  increase  of  the  pore  fluid  pressure  [is  liquefaction]. 
It  involves  a  temporary  transformation  of  the  material  into  a  fluid  mass." 
(A.S.C.E.,  1958)  (American  Geological  Institute,  1972) 

"The  phenomenon  of  the  loss  of  .strength  of  saturated  granular  soils  during 
earthquakes  is  generally  referred  to  as  liquefaction.  The  process  of 
liquefaction  transforms  an  element  of  soil  from  a  state  of  saturated  granular 
solid  to  a  state  of  viscous  fluid."  (Ghaboussi  and  Wilson,  1973) 

"A  phenomenon  in  which  a  cohesionless  soil  loses  strength  during  an 
earthquake  and  acquires  a  degree  of  mobility  sufficient  to  permit 
movements  from  .several  feet  to  several  thousand  feet."  (Seed  and  Idri.ss, 
1971) 

"Complete  liquefaction  —  when  a  soil  exhibits  no  resistance  (or  negligible 
resistance)  to  deformation  over  a  wide  strain  range,  say  a  double  amplitude 
of  20%, 

Partial  liquefaction  —  when  a  soil  first  exhibits  any  degree  of  partial 
liquefaction  during  cyclic  loading. 

Initial  liquefaction  —  when  a  .soil  first  exhibits  any  degree  of  partial 
liquefaction  during  cyclic  loading,"  (Lee  and  Seed,  1967) 

“The  conventional  use  of  the  term  [liquefaction],  as  it  will  be  used 
throughout  the  study,  refers  to  the  phenomenon  which  takes  place  in  a  mass 
of  soil  during  flow  slides.  Liquefaction,  or  flow  failures  of  a  sand  is  caused 
by  a  substantial  reduction  of  its  shear  strength.”  (Castro,  1969) 


As  is  noted,  the  first  five  definitions  refer  to  liquefaction  as  the  loss  of  shear  strength  of  the 
soil  mass.  The  last  definition  by  Castro  refers  to  flow  failures  or  deformations  as  liquefaction, 
which  is  caused  by  soil  losing  its  shear  strength. 

The  controversy  in  literature  results  when  asked  the  question,  how  much  deformation 
constitutes  liquefaction?  Seed  and  Idriss  (1971)  and  Lee  and  Seed  (1967)  indicate  that  any  flow 
deformation  of  the  soil  mass  constitutes  a  liquefaction  condition.  However,  Castro  (1969)  and 
Casegrande  (1976)  feel  that  the  term  “liquefaction”  should  refer  to  a  condition  of  unlimited 
flow  and  that  "cyclic  mobility"  should  refer  to  a  condition  where  the  deformation  is  arrested 
by  a  pore  pressure  reduction.  The  pore  pressure  reduction  is  a  result  of  soil  dilatancy. 


The  differences,  though  only  slight,  do  exist  in  the  literature.  Seed  (1976).  in  his  state-of- 
the-art  address,  notes  these  similarities  and  proposes  the  terms  "initial  liquefaction  with  limited 
strain  potential”  or  "cyclic  liquefaction"  [cyclic  mobility]  to  refer  to  the  condition  where 
limited  strains,  which  may  be  caused  by  liquefaction  of  demser  sands  for  instance,  are  exhibited 
in  a  soil  mass.  This  suggestion  would  hopefully  allow  all  concerned  to  understand  the  basic 
mechanism. 

Further,  Youd  (1975)  also  attempts  to  answer  the  controsersy  with  the  following 
definition: 

“Liquefaction  is  defined  as  the  tramsformation  of  a  granular  material  from  a 
solid  state  into  a  liquefied  state  as  a  consequence  of  increased  pore-water 
pressures.  Solidification  is  defined  as  the  opposite  process,  that  is.  the 
transformation  of  a  granular  material  from  a  liquefied  state  into  a  solid 
state.  Once  liquefied,  a  granular  material  is  free  to  flow,  until  solidification 
occurs.  If  solidification  occurs  after  a  finite  flow  deformation,  the  condition 
is  termed  limited  flow.  If  flow  continues  unbated  under  constant  total  stre.ss, 
the  condition  is  termed  unlimited  flow." 

This  definition  by  Youd  seemingly  accounts  for  both  "liquefaction"  and  "initial 
liquefaction  with  limited  strain  potential"  as  does  Seed’s  (1976). 

DEFINITION  FOR  THIS  STUDY 

The  definition  used  in  this  study  is:  the  transformation  of  a  soil  mass  from  a  solid  state  to 
a  liquefied  state.  This  transformation  is  a  result  of  an  increa.se  in  the  pore  pressure  cau.sed  by 
the  cyclic  loading  of  an  earthquake.  While  soils  at  depth  may  be  said  to  “liciuefy", 
“liquefaction”  is  further  qualified  by  requiring  the  shear  strength  failure  to  reach  the  ground 
surface.  The  term  “liquefaction  potential”  refers  to  a  seismically  induced  pore  pressure  build¬ 
up  at  depth  which  is  capable,  but  not  necessarily  showing  evidence  of  a  lo.ss  of  shear  strength 
at  the  ground  surface.  The  factors  qualifying  “capable"  in  this  definition,  which  are  the 
variables  potentially  allowing  shear  strength  failure  to  reach  the  surface,  will  be  outlined 
further  in  Factors  Affecting  Liquefaction  (Chapter  2). 
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MECHANISM  OF  LIQUEFACTION 

The  basic  cause  of  liquefaction  is  fairly  well  understood  and  accepted.  As  an  earthquake 
occurs  it  creates  waves  of  energs'  (shock  waves)  in  the  bedrock  that  radiate  awa>'  from  the 
.source  in  all  directions.  Seed  and  Peacock  (1971)  suggest  at  a  particular  site,  when  the  bedrock 
is  excited  by  these  waves,  the  .soil  profile  above  the  bedrock  is  set  into  motion  Shear  wases 
propogate  upwards  from  the  bedrock  to  the  ground  surface.  These  shear  wases  cause  eselic 
shearing  stresses  to  deselop  within  the  soil  mass,  which  leads  to  liquefaction.  The  t\|)ieal 
elements  illustrated  in  Figure  2.1(a)  show  the  stresses  applied  during  eartlujuake  loading.  As 
noted  in  part  b  of  the  referenced  figure,  if  a  granular  soil  mass  is  subjected  to  these  shear 
wases,  under  undrained  conditions,  a  build-up  of  excess  hydrostatic  pore  pressure  deselops. 

In  an  effort  to  best  illustrate  the  three  basic  ideas  mentioned  in  the  definitions  of 
liquefaction,  i.e.,  pore  pressure  build-up,  loss  of  .strength,  and  flow  deformation,  it  is  helpful 
to  discuss  Seed's  (1976)  model.  Seed  indicates  that  the  build-up  of  pore  pressure  results  from 
two  interacting  mechanisms.  As  the  cohesionless  soil  mass  is  subjected  to  cyclic  loading,  there  is 
a  tendency  for  the  structure  of  the  soil  to  change  and  to  decrea.se  in  volume.  As  the  .soil  is 
trying  to  compact,  the  load  is  transferred  to  the  pore  fluid  which  results  in  a  build-up  in  pore 
pressure  under  undrained  conditions. 

The  other  mechanism  that  is  associated  with  the  build-up  of  pore  pressure  is  the  soil- 
structure  rebound.  As  the  load  is  transferred  to  the  pore  fluid,  the  structure  exhibits  an  elastic 
rebound  due  to  the  load  release.  The  structure  will  rebound  enough  to  maintain  a  constant 
volume  within  the  soil  system.  The  combination  of  the  volume  decrease,  due  to  cyclic  loading, 
and  the  soil  rebound,  due  to  load  release,  determines  the  amount  of  excess  pressure  that  is 
generated. 

Figure  2,2  is  a  diagram  used  by  Seed  (1976)  to  illustrate  this  point.  Point  A  on  the  void 
ratio  versus  log  pre.s,sure  compre.ssion  curve  is  the  existing  effective  pressure  on  the  soil  clement. 
As  the  .soil  element  is  subjected  to  a  cyclic  shearing  stress,  the  .soil  tends  to  compact  to  poua  3. 
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FIGURE  2.2.  Schematic  lllu.stration  of  the  Mechanism  of  Pore  Pressure  Generation  During 
Cyclic  Loading  (From  Seed,  1976). 
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A.ss(x-iated  with  the  cN'elic  loading  and  soil  compacting  is  the  transfer  of  load  to  the  pore  fluid. 
As  the  load  is  released  from  the  soil  structure,  the  soil  will  tend  to  rebound  along  a  rebound 
curve  to  point  C.  Under  a  drained  condition  the  cyclic  loading  would  cause  an  effect  of 
moving  from  point  A  to  point  B  with  a  net  volume  decrease.  Under  undrained  conditions, 
where  no  volume  change  can  take  place,  the  net  effect  is  moving  from  point  A  to  point  C  with 
a  build-up  of  p(ire  pressure.  The  net  change  in  pore  pressure  is  determined  from  the  initial  and 
final  effective  stres.ses  on  the  soil  element. 

As  described  in  Chapter  1,  Definition  of  Liquefaction,  if  the  build-up  of  pore  pressure 
reaches  the  value  of  the  initial  effective  confining  pressure,  the  soil  will  then  fail  or  liquefy.  At 
this  point  all  resistance  to  deformation  has  been  overcome  and  the  soil  will  deform  under  the 
applied  loads.  The  amount  of  deformation  will  depend  on  the  density-state  of  the  soil  mass.  If 
the  soil  is  in  a  loo.se  condition  then  the  deformation  could  be  unlimited.  However,  if  the  .soil  is 
in  a  medium  to  dense  state  then  the  soil  will  begin  to  dilate.  When  dilation  occurs  the  pore 
pressure  is  reduced  and  continued  deformation  is  arrested. 


FACTORS  AFFECTING  LIQUEFACTION 


GENERAL 

When  determining  the  liquefaction  potential  of  a  soil  profile  there  are  two  aspects  which 
need  examining.  The  first  is  the  magnitude  of  the  shearing  stresses  produced  in  the  profile  by  a 
given  earthquake.  The  magnitude  of  induced  .shearing  stresses  is  a  function  of  the  seismic 
parameters  used  in  the  analysis.  The  second  aspect  is  the  determination  of  the  magnitude  of 
cyclic  shearing  stresses  required  to  initiate  liquefaction  in  the  soil  profile.  This  magnitude  is 
dependent  upon  the  soil  properties  and  other  geotechnical  parameters.  This  portion  of  the 
chapter  will  discu.ss  the  effect  different  parameters  have  on  both  the  induced  and  required 


shearing  .stres.ses. 


INDUCED  SHEARING  STRESSES  —  THE  SEISMIC  PARAMETERS 


General.  The  main  parameters  generated  by  earthquakes  which  are  ^eneralh  discussed  in 
conjunction  with  liquefaction  are: 

•  Intensity  of  ground  shaking 

•  Duration  of  ground  shaking 

•  Magnitude-distance  relationship. 

Before  each  of  these  parameters  are  explained  individualh .  a  brief  discu.ssion  of  their 
relationship  with  each  other  is  required.  The  most  common  combination  of  seismic  pa  ameters 
is  intensity  and  duration  of  ground  shaking.  These  two  factors  are  listed  in  anahsis  procedures 
and  ca.se  studies  by  many  investigators  as  the  principal  seismic  parameters  (Peacock.  19fi8). 
(Seed  and  Idriss,  1971),  (Lee,  1971),  (Christian  and  Swiger.  197.5).  (Seed.  1976).  (Eerritto. 
1977),  (Dobry,  et  al,  1979).  The  third  parameter,  the  magnitude-distance  relationship,  has 
been  used  in  place  of  intensity  and  duration  as  the  basic  seismic  factor  b\  Youd  and  Perkins 
(1977),  and  others. 

Intensity  of  Ground  Shaking.  The  intensity  of  ground  shaking,  or  ground  surface 
acceleration  governs  the  magnitude  of  the  shearing  stresses  that  are  applied  to  the  soil  elements 
(Seed  and  Idri.ss.  1971).  The  higher  the  acceleration  of  the  soil  profile  the  higher  tlu'  shearing 
stresses  that  will  be  induced.  As  illustrated  b\  Peacock  (1968)  in  Figure  2. .3,  under  identical 
confining  pressures,  the  number  of  cycles  reijuired  to  initiate  li(|uefaction  is  decreased  as  the 
peak  shearing  stress  is  increased.  In  fact,  as  will  be  discussed  later,  a  certain  threshold 
acceleration  is  required  to  even  cause  licpiefaction. 

Ground  surface  accelerations  have  been  related  to  magnitude  of  energx  released  and 
di.stance  to  causitive  faults.  Housner  (1964)  produced  a  curve  utilized  for  years,  however,  it 
was  not  universally  applicable  in  that  it  applied  only  to  stiff  soil  sites.  Seed  and  Idriss  (1971) 
have  more  recentl>’  compiled  enough  data  to  create  a  family  of  curves  with  soil  type  and 
stiffness  also  a  variable.  This  relationship  is  demonstrated  in  Figure  2.4.  Housner  and  .Seed 
both  indicate,  however,  that  the  best  method  to  arrive  at  the  intensity  of  ground  acceleration  is 
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FIGURE  2.3.  Effect  of  Confining  Pres.sure  on  Cyclic  Stres.s  Required  to  Cause 
Liquefaction. 


FIGURE  2.4.  Approximate  Relationship.s  Between  Maximum  Acccieration.s  for  Different  Soil 
Conditions  (P'rom  Seed  and  Idriss,  1971). 


to  measure  it  with  instruments  near  a  seismic  source,  not  neces-sarily  from  magnitude-distance 
relationships. 

Becau.se  there  are  only  a  few  recorded  acceleration  histories,  some  studies  have  u.sed 
earthquake  magnitude  and  distance  to  causitive  sources  as  their  parameters  instead  of  the 
e.stimated  ground  acceleration  values.  There  is  much  more  data  on  earthquake  magnitudes  than 
on  earthquake  accelerations,  therefore,  the  larger  data  base  on  earthquake  magnitudes  would 
provide  a  better  solution  for  development  of  empirical  relationships. 


Duration  of  Ground  Shaking.  The  duration  of  intense  ground  shaking  has  a  significant 
effect  on  the  possibility  of  liquefaction.  It  takes  time  for  the  pore  pressure  to  build  to  a 
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magnitude  large  enough  to  overcome  the  confining  pressure  and  reduce  the  soil’s  resistance  to 
deformation.  If  the  duration  of  strong  shaking  is  not  long  enough  then  liquefaction  cannot 
develop.  This  relationship  is  seen  in  Figure  2.3  as,  once  again,  the  number  of  stre.ss  cycles  to 
initiate  liquefaction  is  a  critical  parameter. 

Another  illustrative  example  of  the  importance  shaking  duration  has  on  liquefaction  is 
pointed  out  by  Seed  and  VVilson  (1967).  In  their  analysis  of  the  Turnagain  Heights  land.slide  in 
the  1964  Alaskan  earthquake,  they  report  the  landslide  started  approximately  90  seconds  after 
the  earthquake  began.  These  authors  suggest  that  had  the  earthquake  lasted  onh'  4.5  seconds, 
the  slide  may  never  have  occurred. 

Seed  and  Idriss  (1971)  and  Lee  and  Chan  (1972)  have  attempted  to  standardize  the 
duration  of  ground  shaking  by  utilizing  an  equivalent  number  of  significant  stre.ss  cycles.  It  is 
well  known  that  an  earthcpiake  produces  erratic  stress  cycles  of  varying  frequencies  and 
magnitudes.  These  cycles  are  hard  to  use  in  an  analysis  procedure  becau.se  of  their  non-uniform 
nature.  As  a  result,  the  above  mentioned  investigators  deseloped  a  procedure  where  the  effects 
of  the  significant  earthrjuake  cycles  were  simulated  by  a  certain  number  of  uniform  stre.ss  cycles 
(called  the  ecpiivalent  number  of  significant  stress  cycles).  Each  uniform  stress  cycle  requires  a 
certain  amount  of  time  to  oscillate,  and  so  the  duration  of  strong  ground  shaking  is  determined 
by  the  number  of  uniform  cycles  that  are  applied  to  the  .soil.  The  equivalent  number  of 
significant  stress  cycles  is  used  extensively  in  laboratory  work  and  computer  simulation  to 
model  earthquake  conditions. 

Magnitude-Distance  Relationship.  The  reasons  for  the  development  and  use  of  this 
relationship  have  already  been  discussed  in  the  section  on  Intensity  of  Ground  Shaking.  It 
appears  intuitively  obvious  that  waves  are  attenuated  the  farther  they  travel  and  that  the  shear 
stress  magnitude  subsequently  diminishes  with  the  distance  from  the  causitive  source. 

REQUIRED  SHEARING  STRESSES  -  THE  GEOTECHNICAL  PARAMETER 

General.  The  potential  for  liquefaction  depends,  not  only  on  the  .seismic  activity  of  the 
region  but  also  on  the  condition  within  a  soil  deposit.  Many  soil  deposits  will  not  liquefy. 
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regardless  of  the  inaynitude  aiul  liiiratioii  of  the  eyelie  shearing  stresses  that  are  applied  to  the 
soil.  The  main  faetors  that  are  eotisidered  to  influenee  a  soil  deposit’s  resistance  to  liquefaction 
are  listed  below  . 


•  Soil  Ts'pe 

•  Relative  Densits’ 

•  Initial  Confining  Pressure 

•  Soil  Structure  or  Method  of  Deposition 

•  Age  of  Dep{)sit 

•  Depth  to  Groundwater 

•  Seismic  History 

Some  of  the  governing  factors,  such  as  soil  type,  den.sity,  initial  confining  pressure,  and 
depth  to  groundwater  are  fairly  ssell  understood.  Their  effect  on  the  liquefaction  mechanism 
can  be  recognized  and  quantified.  The  other  factors  have  been  recognized  as  having  some 
effect  on  the  liquefaction  mechanism,  but  they  are  not  fully  understood. 

Soil  Type.  Soil  type  has  been  established  as  one  of  the  main  contributing  factors  to 
liquefaction.  Many  studies  and  investigations  have  pointed  to  the  grain  size  di.stribution  curses 
of  the  soils  involved  as  a  major  factor  in  evaluating  liquefaction  potential.  Fine  uniformly- 
graded  sands  are  cited  as  the  most  susceptible  to  liquefaction.  This  includes,  hossever.  a  range 
from  large  silt  particles  to  medium  coarse  sands  that  could  also  be  classified  as  very  susceptible 
to  liquefaction.  This  range  of  particle  sizes  seems  to  allow  the  build-up  of  excess  pore  pressures 
more  readily  than  any  other  size  di.stribution. 

Coarse  grain  sands  and  gravelly  deposits  have  experienced  .some  cases  of  liquefaction,  but 
in  general  have  a  higher  resi.stance  to  liquefaction  than  the  finer  sands  (Wong,  et  al,  1975), 
(Ross,  et  al.  1969).  The  reason  for  this  is  the  rapid  dissipation  of  exce.ss  pore  pressures.  Gravelly 
deposits  allow  the  excess  pore  pressures  to  dis.sipate  so  rapidly  that  they  do  not  build  up  to  the 
effective  confining  pre.ssures  and  cause  failure. 


Clays,  while  also  being  exempt  from  liquefaction,  do  not  maintain  stal)ilit>  by  pore 
pressure  dissipation  while  loaded.  Their  stability  is  maintained  due  to  their  plastic,  cohesiie 
nature.  Clays  maintain  their  cohesive  properties  with  an  electrostatic  attraction,  rather  than 
the  frictional  intergranular  properties  of  silts  and  sands.  While  friction  is  obviously  reduced  by 
buoyancy,  clays  maintain  their  attraction  through  a  significant  increase  in  excess  pore-water 
pressure.  Seed  (1968)  explains  that  during  earthiiuakes.  the  cyclic  strain  developed  in  elements 
of  sand  in  the  field  mav'  not  be  more  than  .4'';  ,  while  the  strain  in  deposits  of  .soft  clay  mav  be 
as  large  as  1  or  2%.  Figures  2.5  and  2.6  show  the  effects  of  cyclic  strains  of  these  magnitudes 
on  saturated  sand  and  clav'.  In  Figure  2.5  the  application  of  a  cvclic  strain  of  about  .17T  to  a 
sample  of  saturated  sand  caused  initial  li(|uefaction  and  total  failure  after  about  16  strain 
cycles.  In  contrast.  Figure  2.6  shows  the  loss  of  strength  of  a  sample  of  San  Francisco  Bay  mud 
after  200  strain  cycles  (Thiers,  1968).  Even  200  cycles,  with  strain  amplittides  of  up  to  2'~i 
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FIGURE  2.5.  Cyclic  Strain  Test  on  Loose  Sand. 
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FIGURE  2.6.  Effect  of  200  Strain  Cycles  on  the  Strength  of  Clay. 


cau.sed  a  strength  reduction  of  only  about  5%  while  exhibiting  a  continual  deformation.  This 
deformation  is  known  to  continue  as  long  as  the  earthquake  loading  does  not  exceed  the  static 
strength  of  the  clay.  Stress  in  excess  of  this  strength  will  result  in  a  rupture  type  failure. 

Clays,  however,  do  exhibit  another  form  of  dynamic  failure.  This  is  a  result  of  the  collap.se 
of  a  previously  stable  particle  structure  which  is  not  a.ssociated  with  a  pore  pressure  build-up. 
A  clay  susceptible  to  this  phenomenon  is  termed  sensitive  and  thixotropic. 

A  new  development,  even  for  the  most  liquefiable  soils,  are  the  grain  characteristics. 
Several  studies  (Annakai.  1975),  (Castro  and  Poulos.  1976)  have  shown  that  different  types  of 
sands  with  essentially  the  same  grain  size  curves,  compacted  to  the  same  relative  density,  and 
compacted  with  the  same  compaction  methods,  differ  as  to  their  liquefaction  potential.  This 
phenomenon  is  not  completely  understood  as  to  what  causes  the  observed  difference,  but  is 
probably  related  to  grain  size. 

Relative  Density.  Relative  density  has  long  Ireen  recognized  as  a  major  factor  affecting  the 
liquefaction  potential  of  a  deposit.  If  a  deposit  is  in  a  relatively  loose  condition,  i.e..  low 
relative  density,  licjuefaction  can  be  initiated  by  lower  .shearing  stresses  or  shorter  durations  of 


shaking  than  if  the  deposit  is  at  a  higher  relative  density.  Figure  2.7  illustrates  results  of  a 
typical  cyclic  load  test  run  on  loose  sand  (Seed  and  Peacock,  1968).  In  this  test  a  cyclic  shear 
stress  of  constant  amplitude  (  ±  .33  kg.  per  ,sq.  cm.)  was  applied  with  a  frequency  of  one  cycle 
per  second  under  a  vertical  con.solidation  load  of  5  kg.  per  sq.  cm.  The  resulting  changes  in 
horizontal  shear  strain  and  pore- water  pressure  were  recorded. 
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FIGURE  2.7.  Pulsating  Load  Test  on  Loose  Sand — Simple  Shear  Conditions. 


It  is  interesting  to  note  that  in  the  initial  stages  of  the  test  (i.e.,  the  first  24  cycles  of  stress 
application),  the  gradual  increase  in  pore- water  pressure  did  not  produce  a  significant  increase 
in  the  shear  deformation,  even  though  the  effective  stress  was  con.stantly  being  reduced. 
However,  as  soon  as  the  effective  stress  within  the  sample  had  been  reduced  to  zero,  thereby 
reducing  the  shear  strength  to  zero,  the  deformations  increased  rapidly  and  failure  due  to 
liquefaction  occurred  within  2  cycles.  This  sudden  failure  appears  to  be  characteristic  of  loose 
sands  in  a  variety  of  dynamic  te.sts. 
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In  contrast.  Figure  2,8  exhibits  results  of  the  same  cyclic  load  test  on  dense  sand  (Seed  and 
Peacock,  1968).  An  analysis  of  the  pore-uater  pressure  and  shear  strain  curses  shoss'  that  dense 
sands  exhibit  somewhat  the  same  bchasior  as  loo.se  sands.  Initially,  there  was  no  noticeable 
deformation  with  increasing  pore-ssater  pre.ssure.  Hosseser.  unlike  loose  sand,  the  strain 
amplitude  increased  more  slowly  with  increasing  number  of  stre.ss  cycles.  Thus,  the  dense  sand 
exhibits  a  gradual  build-up  of  shear  strain  with  the  onset  of  initial  liquefaction,  but  shear 
strains  tend  toward  a  limiting  salue  of  about  8‘t  double  strain  amplitude.  A  conclusion 
tspified  in  numerous  tests  is  that  dense  sand  is  capable  of  undergoing  "initial  liquefaction  with 
a  limited  strain  potential”  (Seed,  1976) 


FIGURE  2,8.  Pulsating  Load  Test  on  Dense  .Sand — Simple  Shear  Conditions. 


The  theories  for  such  a  drastic  difference  in  the  behavior  of  sands  and  .silts  of  different 
relative  density  is  still  being  debated.  As  pointed  out  by  Norris  (1980),  the  arguments  between 
Casegrande  and  Seed  are  growing  closer  together.  Norris  .suggests  that  Casegrande’s  critical 
void  ratio  theory  seems  to  be  a  good  starting  point  for  debate.  As  discussed  in  Mechanisms  of 
Liquefaction,  the  critical  void  ratio  suggests  a  void  ratio  where  values  le.ss  than  this  ratio  (i.e.. 
in  a  denser  state)  tend  to  dilate  during  shear  strain  as  the  particles  are  forced  to  "ride  up  and 
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()\er  each  other.  ”  In  contrast,  void  ratios  in  e.xcess  of  this  \alue  (i.e..  in  a  looser  state)  tend  to 
collapse  into  a  denser  state  during  shear  .strain  thereby  tending  toward  a  total  solnrne  decrease. 
Clearly  a  decrease  in  total  soliime  \sill  result  in  an  increase  in  pore-water  pressure  when 
sample  drainage  and  subsecpient  actual  volume  change  is  prevented. 

In  dense  sand  and  silt,  however,  as  the  effective  normal  stress  approaches  zero  and  a 
significant  amount  of  shear  strain  is  initiated,  the  tendency  is  to  cause  a  volume  increa.se, 
therebs'  causing  a  tensional  pore-water  pressure  build-up  in  the  undrained  condition.  As  the 
shear  strain  proceeds,  the  e.xcess  pore-water  pre.ssure  builds  in  a  negatise  direction  and  the 
normal  effective  stress  increa.ses  (o'  =  a  -  (-ft  excess)).  As  a  result,  with  each  cvcle  only  a 
certain  amount  of  deformation  can  occur  before  sufficient  resistance  is  mobilized,  therebs' 
realizing  Seed's  “limited  strain  potential.” 

It  is  with  regard  to  the  effective  stress  approaching  zero,  and  whether  dense  sands  actually 
liqueh'  that  Casegrande  and  Seed  founded  their  debate,  Casegrande,  using  his  critical  void 
ratio  concept,  argued  that  an>  sand  existing  at  a  void  ratio  below  the  critical  \alue  tends  to 
dilate  during  dynamic  shear  stress,  thereby  negating  any  pore  pressure  build-up  as  the  sample 
would  be  in  tension.  Seed,  on  the  other  hand,  argues  that  the  critical  void  ratio  is  applicable  to 
static  loading,  but  that  it  does  not  deal  with  the  po,s,sible  volume  decrease  tendency  of  dense 
sand  during  the  pore  pressure  build-up  pha.se  occurring  prior  to  the  shear  strains  causing 
dilation. 

Initial  Confining  Pressure.  The  confining  pressure  that  a  soil  element  is  under  is  of 
significant  importance  in  determining  its  liquefaction  potential.  The  higher  the  effective 
confining  pre.ssures  the  higher  the  excess  pore  pressures  need  to  be  to  overcome  the  strength  of 
the  soil  that  prevents  deformation. 

Regarding  the  overburden  pressure,  the  deeper  the  liquefiable  material  is  in  the  profile, 
the  higher  the  overburden  pressure  and  the  more  resistant  the  .soil  becomes  (Seed  and  Lee, 
1966).  Meslov  (1957)  and  Ambra.seys  and  Sarma  (1969)  u.se  this  to  suggest  that  the  presence  of 
a  surcharge  will  reduce  the  tendency  of  a  deposit  to  liquefy. 
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Anothfr  itnportaiil  cdiifinint'  pressure,  the  lateral  earth  pressures,  can  he  a  function  of 
either  present  overburden  pressures  or  the  geologic  history  of  the  deposit.  If  a  deposit  has  been 
svibjected  to  higher  overburden  pressures  than  now  e.xists  on  the  deposit,  it  has  been 
overconsolidated  to  a  degree,  B\-  increasing  the  overconsolidation  ratio  (ratio  of  highest  past 
overburden  pre.ssure  to  present  overburden  pressure),  an  increase  in  lateral  earth  pressure  is 
also  produced.  The  result  of  increasing  the  lateral  pressure  is  the  same  as  increasing  the 
overburden  pressure,  a  more  resistant  soil  deposit  to  liquefaction  (Seed  and  Peacock,  1971). 

Lemuk  (1966)  and  Shannon  (1966)  illustrate  the  effects  of  initial  confining  pressures  as 
they  de,scribe  the  Seward,  Alaska  flowslide  utilizing  eyewitness  accounts  coupled  with  a 
detailed  geotechnical  analysis.  In  short,  liquefaction  occurred  at  depths  of  50  feet  or  less 
causing  flowslides  downhill  into  the  harbor.  When  one  such  layer  slid  away,  the  reduction  in 
confining  pressure  permitted  liquefaction  of  the  underlying  layer  and  this  procedure  could 
repeat  progre.ssively  during  the  earthquake. 

Seed,  Lee,  and  Idriss  (1969)  suggest  their  laboratory  tests  show  that  the  larger  the  ratio  of 
initial  shear  stress  to  initial  normal  confining  pressures  acting  on  a  horizontal  surface  of  a  soil 
element,  the  greater  is  the  horizontal  cyclic  .shear  .stress  required  to  induce  liquefaction  in  a 
given  number  o*^  stress  cycles.  This  relationship  leads  directly  to  the  conclusion  that  other 
factors  being  equal,  liquefaction  will  be  induced  more  easily  under  le\el  ground  conditions 
than  is  sloping  zones  of  a  deposit. 

The  parameter  that  many  laboratory  studies  base  their  results  on  is  the  cyclic  stre.ss  ratio. 
The  cyclic  stress  ratio  is  the  shearing  .stre.ss  required  to  cau.se  liquefaction  to  the  effective 
confining  stress  on  the  sample.  The  cyclic  stress  ratio  is  usually  plotted  versus  the  number  of 
equivalent  uniform  stre-ss  cycles  that  is  required  to  cau.se  liquefaction.  The  shearing  stress  is 
proportionate  to  the  confining  stress  so  the  cyclic  stress  ratio  provides  a  consenient 
dimensionless  parameter  that  includes  two  factors  that  influence  the  liquefaction  characteristics 
of  a  .soil  deposit. 


Soil  Structure  (Method  of  Deposition).  Soil  structure  and  the  effect  that  it  has  on  the 
liquefaction  characteristics  of  a  .soil  is  a  fairly  recent  finding.  Mulilius.  et  al  (1975)  clearU'  show 
that  the  method  in  which  a  soil  depo.sit  has  been  laid  down  makes  a  difference  in  its 
li<luefaction  characteristics.  Their  report  deals  mainly  with  different  methods  of  preparing 
laboratory  samples,  but  the  conclu.sions  are  easily  extrapolated  to  field  conditions.  This  means 
that  a  cohesionless  depo.sit  probabb'  will  have  different  potentials  for  litjuefaction  depending  on 
whether  it  was  deposited  b\'  flusial  deposition,  direct  sedimentation,  or  by  eolian  deposition. 
Each  soil  structure  would  be  different  and  would  produce  different  susceptibilities  to 
liciuefaction. 

The  exact  nature  of  what  the  soil  structure  does  to  alter  litjuefaction  potential  is  not 
completeh  understood.  It  is.  therefore,  hard  to  (juantify  and  to  indicate  hou  the  soil  structure 
could  be  considered  in  an  analysis  procedure. 

Age  of  Deposit.  The  length  of  time  that  a  soil  depo.sit  has  been  in  place  also  has  an  effect 
on  its  liciuefaction  characteristics.  The  older  a  soil  deposit  is  the  less  chance  there  seems  to  be 
that  it  will  licpiefy.  Holocene  (recent)  and  late  Pleistocene  deposits  are  cited  as  the  most  likely 
deposits  to  liquefs'.  This  conclusion  is  based  on  numerous  ca.se  studies  that  determine  the  age  of 
deposits  that  were  known  to  li<iuef\  during  a  seismic  event  (Youd  and  House.  1977).  (Youd  and 
Perkins,  1977).  Recent  laborators  studies  have  also  shown  that  the  longer  a  sample  is  allowed 
to  sit  before  testing  the  more  resistant  the  samples  becomes  to  liquefaction.  Lee  (1975)  indicates 
that  this  increa.se  in  resistance  might  be  the  result  of  cementation  between  the  contact  points  of 
sand  grains,  or  the  development  of  a  more  stable  structure  resulting  from  secondars' 
compression. 

Depth  to  Groundwater.  The  depth  to  the  groundwater  level  plays  a  major  role  in 
liquefaction  susceptibility.  If  a  soil  is  not  fully  saturated  then  it  is  impossible  to  develoii  excess 
pore-water  pressure. 

The  depth  to  the  water  table  also  affects  the  confining  pressure  on  the  soil  element.  The 
higher  the  water  level  in  the  soil  profile,  the  lower  will  be  the  effective  confining  pressure  at 
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any  depth  below  the  water  le\el.  This  indicates  that  a  high  water  table  in  a  soil  deposit  not 
only  saturates  the  deposit,  making  licpiefaction  possible,  but  also  reduces  the  effectise  pressures 
on  the  soil  elements  below  the  water  le\el.  If  the  effective  confining  pressures  are  reduced  on 
the  soil  elements,  then  the  deposit  is  more  su.sceptible  to  lif|uefaction. 

Seismic  History.  Although  not  comirletely  understood,  seismic  or  strain  hi.story  has  a 
significant  effect  on  the  susceptibility  of  a  soil  deposit  to  li(}uefaction.  Seed,  et  al  (1977)  and 
Youd  (1977)  states  that  this  change  in  liquefaction  susceptibility  could  possibly  result  from  a 
volume  change,  i.e.,  interparticle  packing  or  structural  changes  caused  by  previous  earthquakes 
that  change  the  pore  [iressure  build-up  mechanism  or  ability.  The  effects  were  first  pre,sented 
by  Finn,  et  al  (1970).  but  later  studies  have  shown  some  of  the  same  results  (Lee  and  Focht, 
197.5).  (Seed,  et  al,  1977),  (Youd.  1977). 

Youd  (1977)  summarizes  the  results  by  .stating  4  conclusions.  (1)  Drained  prestraining 
generally  produces  a  packing  in  sandy  materials  that  is  more  resistant  to  li(juefaction  than  was 
the  original  packing.  (2)  Undrained  prestraining,  not  producing  liquefaction,  with  subsecjuent 
pore-pressure  rvdief  also  produces  a  packing  that  is  more  resistant  to  liquefaction  than  the 
original  packing.  (3)  Undrained  prestraining  producing  liquefaction  and  limited  flow 
deformation  leaves  a  packing  highly  susceptible  to  liquefaction  even  after  recon.solidation.  (4) 
Following  a  loading  phase,  which  produces  liriuefaction,  and  a  limited  flow  deformation  and 
reconsolidation  phase  as  outlined  in  number  3,  additional  drained  loading  cycles  produce 
significant  increases  in  the  liquefaction  resistance  of  previously  liquefied  sand  (Y'oud,  1977). 
Suggested  causes  of  increased  liquefaction  susceptibility  with  previously  liquefied  soil  include  a 
loss  of  strengthening  from  cementation,  prior  strain  history,  and  lost  grain  structure. 

METHOD  OF  ANALYSIS 


GENERAL 

With  all  the  information  now  being  generated  concerning  liquefaction,  some  authors  have 
attempted  to  form  a  rational  method  of  analy.sis.  Each  different  method  is  ba.sed  on  the  seismic 
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and  geotechnical  parameters  previoush  discussed.  Individual  authors  mat’  weigh  individual 
parameters  differently,  depending  on  their  feelings  of  importance  to  the  situation.  \o  one 
method  is  clearl.v  better,  nor  has  one  method  been  universall.v  accepted.  Onh’  with  numerous 
trial  and  error  situations,  with  the  final  test  being  actual  case  histories,  will  the  most  accurate 
method  of  analysis  be  accepted.  Currently,  some  methods  might  be  used  more  than  others,  but 
it's  generally  only  because  thev  have  been  in  existance  longer, 

LABORATORY  TESTING  PARAMETERS 

A  large  amount  of  laboratorv'  work  has  been  carried  out  by  investigators  in  the  analysis  of 
the  dvnamic  stress  field  in  .soil  during  liquefaction.  .Manv'  of  the  testing  procedures,  together 
with  the  subse()uent  results,  have  been  used  as  a  springboard  for  other  field  methods  for 
analyzing  liquefaction  potential.  This  laborator>  testing  section  is  meant  to  (trovide  an 
understanding  of  the  stre.ss  field  and  discuss  two  of  the  more  major  soil  dvnamic  tests  and  the 
potential  problems  of  each.  It  must  first  of  all  be  realized  that  whenever  laboratorv  testing  is 
carried  out  on  soil,  disturbances  of  the  in-situ  soil  structure  is  inevitable  and  therefore 
understood  to  be  a  problem. 

It  has  previously  been  suggested  (Seed  and  Peacock,  1971)  that  in  manv-  ca.ses  it  appears 
that  the  main  forces  acting  on  soil  elements  in  the  field  during  earthquakes  are  those  resulting 
from  an  upward  propagation  of  shear  motions  from  underlying  rock  formations.  The  result  is 
that  a  soil  element  in  a  deposit,  with  a  horizontal  grtjund  surface,  is  subjected  to  the  sequence 
of  stress  application  illustrated  in  Figure  2.9.  Before  the  earthquake  the  soil  element  is 
subjected  to  the  principal  stress  (o')  and  kf,a',  in  which  kp  is  the  coefficient  of  lateral  earth 
pressure  at  rest.  During  the  earthquake  the  element  will  be  deformed  through  the  development 
of  a  .shear  stress  (rj^y)  which  will  reverse  directions  many  times  during  the  cyclic  loading. 

Utilizing  analytical  methods,  such  as  the  program  SHAKE  (Schnabel,  et  al,  1972)  which 
predicts  ground  motion  acceleration,  the  values  of  rj,y  and  a'  can  be  evaluated  for  different 
soil  elements  in  the  profile.  The  magnitude  of  this  shear  stress  will  vary  .somewhat  erratically 
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FIGURE  2,9.  Cyclic  Stress  Loading — Simple  Shear  Test. 


from  one  cycle  to  the  next.  As  mentioned  in  the  discussion  on  Duration  of  Shaking,  by 
appropriately  weighing  the  different  cyclic  stress  amplitudes  this  loading  pattern  can  be 
converted  to  an  equivalent  number  of  uniform  stress  cycles.  It  is  this  pattern  to  which  the 
samples  are  subjected  during  laboratory  testing.  The  objective,  therefore,  of  the  test  procedure 
is  to  determine  the  magnitude  of  the  cyclic  shear  stress  (r|,v)  '"’hich  will  cause  liquefaction  of 
representative  samples  for  different  values  of  the  initial  effective  normal  stress  (a^. ). 

CYCLIC  SIMPLE  SHEAR  TEST 

Seed  and  Peacock  (1971)  note  that  the  laboratory  cyclic  simple  shear  test  is  considered  to 
correctly  simulate  the  field  loading  condition  (Figure  2.9).  Both  field  and  laboratory  elements 
are  first  coasolidated  under  an  effective  overburden  pressure  (a(. ).  and  because  the  sample  is 
restrained  (see  Figure  2.10  for  test  apparatus)  from  lateral  deformation,  a  lateral  earth  pressure 
is  modeled  equal  to  kgOy ,  yielding  the  stress  condition  as  shown  in  Figure  2.9.  During  an 
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FIGURE  2.10.  Laboratory  Cyclic  Simple  Shear  Test  Apparatus. 


earthquake,  or  as  a  result  of  stress  applications  in  the  laboratory,  the  soil  element  is  subjected 
to  a  horizontal  and  vertical  shear  stress  (r|,y).  This  stress  condition  changes  as  shown  by 
Figure  2.9’s  accompanying  Mohr  s  circle.  Note  that  as  T},y  is  applied  in  a  csclic  manner,  the 
maximum  shear  stress  (r^^x)  sample  increases,  and  in  fact  is  always  greater  than  rj,y 

when  kg  is  less  than  1,  Note  also  that  the  .shear  stress  required  to  cau.se  failure  increases  with 
increasing  k^,  or  confining  pressure.  This  further  illustrates  the  importance  of  confining  stress. 

It  has  therefore  been  established  that  the  application  of  stress  in  the  cyclic  simple  shear 
test  is  similar  to  that  in  the  field.  Note,  however,  that  when  modeling  liquefaction  phenomena 
in  the  laboratory,  testers  are  interested  in  working  at  or  near  a  failing  eondition  at  all  times. 
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From  this  we  can  clearly  see  that  the  effects  of  nominiforni  hehavinr  can  he  \  erv  detrimental  to 
the  results. 

Seed  and  Peacock  point  out  that  accurate  simulation  of  field  stress  conditions  are  limited 
by  the  following  factors: 

1.  Because  of  the  difficulties  in  enclosing  test  .specimens  to  maintain  constant 
s'olume  or  uniform  strain,  it  is  difficult  to  prepare  test  specimens  in  a 
uniform  condition,  i.e.,  uniform  specimen  densities,  representatise  of  field 
conditions. 

2.  Difficulties  are  encountered  in  ensuring  a  uniform  application  of  shear 
stress  across  the  width  of  the  sample,  the  uniformity  depending  to  a 
considerable  e.xtent  on  the  stress  transfer  mechanism  between  the  cap.  base, 
and  the  sample. 

3.  Roscoe  (1953),  who  has  worked  to  refine  the  testing  ap|)aratus.  notes 
difficulties  in  developing  complementary  shear  stresses  along  the  vertical 
faces  of  a  test  specimen;  the  absence  of  which  k'ads  to  nommiformities  in 
stress  and  .strain  conditions  near  the  edges  of  the  sam\)le  and  to  some 
rocking  of  the  cap  and  base  under  cyclic  shear  stress  application. 

4.  Unle.ss  the  specimen  is  specially  constrained,  concentrated  zones  of  stress 
and  subsequent  strain  are  likely  to  develop  at  the  points  of  a]splication,  i.e.. 
the  cap  and  base. 

Seed  and  Peacock  have  found  that  the  net  effect  of  the  limitations  on  both  the  apparatus  and 
testing  procedures,  and  the  subsequent  nonuniformities  in  .stress  conditions  cause  samples  to  fail 
up  to  35 ‘t  lower  than  would  be  recjuired  in  the  field. 


CYCLIC  TRIAXIAL  TEST 

Because  of  its  more  widespread  availability  and  simplified  procedures,  the  cyclic  triaxial 
compression  test  is  a  popular  tool  for  use  in  evaluating  liciuefaction  potential.  It  must  be  noted, 
however,  that  the  stress  conditions  in  this  test  are  not  readily  comparable  to  those  for  normal 
earthquake  loading  conditions  (Seed  and  Peacock,  1971).  This  difference  in  effective  stresses 
can  best  be  exemplified  by  looking  at  the  variations  in  total  stresses  on  a  particular  plane  at 
inclination  a.  Figure  2.11(a)  illustrates  the  method  of  loading  employed  by  the  cyclic  triaxial 
test  with  a  superimposed  plane  at  inclination  a.  The  accompanying  Mohr  circles  (Figure 
2.11(b))  show  the  stress  directions  involved.  In  comparing  times  1  and  2,  as  also  seen  in  actual 
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FIGURE  2.11.  Cyclic  Stress  Loading — Cyclic  Triaxial  Test. 


field  conditions,  the  shear  stress  magnitudes  are  equal  but  their  directions  are  opposite.  Unlike 
the  actual  field  conditions  we  wish  to  model,  note  the  variations  in  total  normal  stress 
(K'curring  between  the  two  intervals.  As  noted  in  the  discussion  of  the  cyclic  simple  shear  test,  a 
typical  seismically  loaded  field  element  is  subjected  to  a  relatively  constant  normal  load  with 
the  cyclic  shear  component  being  opposite  in  direction  but  equal  in  peak  magnitude.  Figure 
2.11  suggests  that  only  when  a  =  45°  (or  2a  =  90°)  will  the  modeled  laboratory  stress 
conditions  be  similar  to  those  expected  in  the  field. 

Seed  and  Peacock  (1971)  caution  that  it  is  only  possible  to  produce  the  desired  stress 
changes  isotropically  by  con.solidating  the  sample,  i.e.,  under  a  lateral  pressure  coefficient 
equal  to  1  (k^  =  1).  If  any  other  initial  stre.ss  condition  is  used,  there  is  no  plane  in  the  sample 
which  receives  the  desired  symmetrical  changes  in  shear  stress  corresponding  to  those  on  the 
horizontal  yiane  in  the  field.  This  factor  further  alienates  the  triaxial  model  from  the  desired 
field  conditions  as  the  accepted  in-situ  k^,  value  usually  ranges  from  .4  to  .5. 
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Seed  and  Peacock's  (1971)  further  observations  regarding  the  cyclic  triaxial  test  include  the 
following  factors  potentially  affecting  the  accuracy; 


1.  Under  field  loading  conditions,  the  direction  of  the  major  principal  stress 
will  \ar\  from  about  0“^  to  40°  on  each  side  of  the  vertical,  but  in  the 
triaxial  compression  test  the  direction  of  the  major  principal  stress  max- 
rotate  through  an  angle  of  90  on  each  stress  cycle. 

2.  Unless  special  precautions  are  taken,  friction  between  the  cap  and  base  of 
the  test  specimen  can  cau.se  stress  concentrations  which  can  lead  to 
P’-emature  failure  of  the  specimen  (Castro,  1969). 

3.  During  the  lateral  compression  part  of  the  stress  cycle,  the  sample  will 
tend  to  "neck"  as  failure  approaches,  and  it  is  difficult  to  determine  the 
stress  condition  on  the  sample  once  this  occurs  (i.e..  the  area  is  constantlx’ 
changing).  The  effect  is  not  likely  to  be  very  .significant  for  relativelx  loo.se 
samples  in  which  licpiefaction  and  the  development  of  large  strains  occur 
almost  simultaneously,  but  it  can  have  a  significant  effect  on  the  behavior 
of  dense  samples,  which  in  turn  would  lead  to  apparentlv-  larger 
deformation  than  would  otherwi,se  be  the  case. 

4.  There  appears  to  be  evidence  that  .stresses  causing  li(iuefaction  in  triaxial 
compression  tests  are  influenced  by  the  diameter  of  the  specimen,  where 
the  smaller  diameter  specimens  require  higher  .stres.ses  to  cause 
liciuefaction. 

The  effects  of  the  first  three  limitations  are  likely  to  cau.se  test  specimens  to  licjuefv  or  fail 
at  stresses  lower  than  those  which  would  be  recjuired  under  corresponding  field  conditions. 
However,  Seed  and  Peacock  (1971)  al.so  advise  that  utilizing  a  k^  =  1  condition  in  the 
laboratory  will  cause  stres.ses  inducing  failure  or  liquefaction  to  be  higher  than  those  inducing 
failure  in  the  field. 


As  an  increasing  amount  of  laboratory  analysis  is  undertaken  regarding  liquefaction,  the 
relationships  between  cyclic  simple  shear  tests,  cyclic  triaxial  tests,  and  actual  field  conditions 
are  constantly  being  updated.  It  is  recognized  that  there  are  numerous  problems  in  attempting 
an  accurate  field  model.  The  employment  of  correction  factors,  the  derivations  of  which  are 
beyond  the  scope  of  this  paper,  are  now  being  received  w-ith  an  increasing  amount  of 
confidence  as  the  potential  for  catastrophes  due  to  liquefaction  is  being  realized. 


.39 


SIMPLIFIED  PROCEDURE 

Over  a  period  of  time,  the  experienced  gained  from  the  rigoroas  anaKsis  employing  the 
shear  stress  histories  computed  via  SHAKE  (Schnabel,  ct  al,  1972),  or  any  other  appropriate 
program  has  led  Seed  and  Idriss  (1971)  to  propose  a  simplified  method  which  eliminates: 

1.  the  need  to  determine  the  variation  in  shear  stres.s  history  with  depth,  and 

2.  the  need  to  perform  laboratorv  te.sts. 

This  method  is,  of  course,  not  as  accurate  as  the  more  rigorous  analysis.  If  indications 
from  the  method  reveal  no  danger  of  liciuefaction  by  a  wide  margin,  a  sate  condition  can  be 
indicated.  If  the  procedure  is  either  safe  by  a  narrow  margin,  or  unsafe,  a  more  detailed 
analysis  should  be  employed. 

This  method  is  operational  due  to  some  simplifying  assumptions.  The  basic  premise  is  to 
compare  the  cxclic  stres.ses  that  an  earthquake  will  cause  in  a  soil  profile  to  the  cyclic  stres.ses 
that  are  required  to  cau.se  liquefaction  in  that  same  profile,  as  illustrated  in  Figure  2.12.  The 
overlap  region  on  the  figure  is  the  area  of  conc-ern. 

The  stresses  induced  b\'  an  earthquake  are  estimated  by  a  simplified  equation  of  the  form; 

W  =  -65  Y  '■d  (2.1) 

where 

Tave  =  average  shearing  stress  caused  by  the  earthquake 

7  =  unit  weight  of  the  .soil 

g  =  acceleration  of  gravity 

^max  ~  maximum  ground  surface  acceleration 

rd  =  stress  reduction  coefficient 

The  average  shearing  stress  is  based  on  the  amount  of  stress  that  will  be  realized  beneath  a 
rigid  column  of  .soil  at  a  depth  h.  The  stress  reduction  coefficient,  rd,  is  used  becau.se  the  soil 
column  is  not  truly  rigid.  The  multiplier  constant  is  an  assumption  that  the  average  stress  is 
65  %  of  the  maximum  stress  induced  by  the  irregular  stress  history.  Figure  2.12  is  based  on 
numerous  calculations  of  equivalent  uniform  .shearing  stresses  for  different  stress  histories.  The 
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FIGURE  2.12.  Method  of  Evaluating  Liquefaction  Potential  (From  Seed  and 
Idriss,  1971). 


maximum  ground  surface  acceleration  is  taken  from  a  family  of  curves  proposed  by  Seed,  who 
integrated  a  known  bedrock  acceleration  with  a  known  type  of  soil  profile  to  yield  a  surface 
acceleration.  (See  Figure  2.4) 
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The  duration  of  ground  shaking  is  accounted  for  by  adjusting  the  number  of  equivalent 
uniform  cycles  that  are  assumed  to  be  applied  to  the  soil  profile  by  the  earthquake.  This  value 
can  be  obtained  from  Seed  and  Idriss'  recommendation  which  is  totalK'  a  function  of 
earthquake  magnitude. 

Shearing  stre.sses  retjuired  to  cause  liquefaction  in  the  soil  profile  are  usually  determined 
on  the  basis  of  laboratois'  tests.  Dynamic  triaxial  shearing  tests,  like  those  previously  described, 
are  usually  run  to  determine  the  cyclic  stress  ratio  required  to  cau.se  liquefaction  in  a  gi\en 
number  of  cycles  and  at  a  given  relative  density.  Stress  ratio  also  depends  on  the  mean  grain 
size  diameter,  D5Q. 

If  dynamic  tests  are  not  performed,  then  a  cyclic  stress  ratio  can  be  estimated  for  a  gi\en 
D5Q  value  at  a  given  relative  density  and  at  a  given  number  of  uniform  stress  cscles  from 
dynamic  test  data  run  on  other  samples.  The  form  of  the  equation  used  for  estimating  stresses 
required  to  cause  liquefaction  is 


where 


(2.2) 


1  Dr  =  cyclic  stress  ratio  causing  liquefaction  at  a  relative  density  of  Dr 


ratio  of  the  deviator  stress  to  the  initial  ambient  pressure  that  cau.ses 
liquefaction  at  a  relative  density  equal  to  Dr';  from  cyclic  triaxial  test. 


Cr 

Dr 

Dr' 


correction  factor  to  correct  triaxial  data  to  field  conditions. 

relative  density  ratio,  to  change  data  from  a  relative  density  of  Dr'  to 
a  given  relative  density  of  Dr  (Dr'  =  50%  for  Seed’s  data). 


This  form  of  the  equation  is  valid  for  relative  densities  to  80%  . 


If  the  average  .shearing  stress  (I'ave)  from  equation  2.1  is  set  equal  to  the  shearing  stress  (r) 
from  equation  2.2,  then  the  maximum  acceleration  r'an  be  determined.  Both  equations 

must  relate  to  the  same  number  of  shearing  stress  cycles.  The  value  of  can  be  plotted 

versus  the  value  of  Dr  for  that  .set  of  data.  Different  relative  densities  can  be  used  to  arrive  at 


I 
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different  values  of  plot  forms  a  boundary  between  liquefaction  and  non-liquefaction 

conditions.  An  example  of  this  type  of  plot  is  shown  in  Figure  2.13, 


FIGURE  2.13.  Evaluation  of  Liquefaction  Potential  for  Very  Fine  Sand — 20  Stress 
Cycles  (From  Seed  and  Idriss,  1971). 


EMPIRICAL  METHODS 

Another  method  that  is  becoming  widely  known  is  presented  in  its  most  recent  form  by 
Seed,  et  al  (1977).  This  method  is  based  on  empirical  data  of  .sites  that  have  been  studied 
where  liquefaction  has  or  has  not  occurred.  The  cyclic  stress  ratio  causing  liquefaction  was 
plotted  versus  blow  count  data  from  the  standard  penetration  test,  corrected  to  an  effective 
overburden  pressure  of  1  ton/sq.  ft.  A  lower  bound  curve  was  established  that  separates  the 
liquefaction  conditions  from  the  non-liquefaction  conditions,  as  depicted  in  Figure  2.14. 


CYCLIC  STRESS  RATIO  CAUSING  LIQUEFACTION,  Til), 
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FIGURE  2.14.  Correlation  Between  Stres.s  Ratio  CaiLsing  Liquefaction  in  the  Field  and 
Penetration  Resi.stance  of  Sand  (From  Seed,  Mori,  and  Chan,  1977). 
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To  use  the  figure,  the  cyclic  stress  ratio  must  be  determined.  First,  the  shearing  stresses 
that  would  be  created  in  a  soil  profile  by  an  earthquake  need  to  be  predicted.  This  can  be 
accomplished  using  the  same  equation  2.1  as  given  for  the  simplified  procedure  for  evaluating 
soil  liquefaction  potential.  Then,  the  initial  effective  confining  stress  used  to  form  the  ratio  can 
be  determined  from  boring  data  taken  at  the  site.  The  ratio  of  the  shearing  stress  to  the  initial 
confining  stress  is  plotted  on  the  chart  versus  the  corrected  blow  count  also  obtained  from 
boring  data.  If  the  point  is  above  the  boundary  line  then  liquefaction  is  a  possibility.  If  the 
point  plots  below  the  line  then  liquefaction  probably  will  not  develop.  If  the  point  plots  close 
to  the  line,  on  either  side,  then  a  closer  look  at  the  data  using  a  different,  more  detailed 
method  of  analysis  is  probably  justified.  This  type  of  approach  could  be  used  to  pick  out  the 
sites  and  profile  layers  that  might  need  more  attention. 

Whitman  (1971)  presented  a  method  very  similar  to  the  above  empirical  method.  He  also 
used  data  from  earthquakes  that  have  caused  liquefaction,  as  well  as  from  a  few  that  did  not. 
The  data  that  he  plotted  was  the  stress  ratio  versus  relative  densits .  He  pointed  out  that  the 
data  is  not  sufficient  to  define  a  trend  or  a  boundary  line,  but  w’ith  more  information  on  other 
earthquakes  some  type  of  distinction  could  be  made  between  liquefiable  and  non-liquefiable 
deposits. 

The  main  difference  between  Whitman’s  chart  and  Figure  2.14  is  the  plotting  of  relative 
densities  instead  of  blow  count  data.  The  relative  densities  are  determined  using  standard 
penetration  results  and  relationships  from  Gibbs  and  Holtz  (1957).  It  is,  however,  more 
appropriate  to  plot  the  standard  penetration  results  rather  than  the  relative  densities.  The  blow 
count  data  in  some  ways  accounts  for  more  of  the  factors  that  influence  a  soil's  liquefaction 
potential  than  does  relative  density  alone.  As  the  factors  such  as  relative  density,  .soil  particle 
cementation,  lateral  earth  pressure,  and  prior  .seismic  hi.stories  affect  liquefaction  potential  they 
similarly  affect  blow  count  data.  It  is,  therefore,  suggested  that  the  blow  count  data  gives  a 
better  representation  of  true  strength  resulting  from  many  influencing  factors  and  not  just  one 
(Seed,  et  al,  1977). 
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Donovan  and  Singh  (1976)  presented  a  method  for  determining  liquefaction  susceptibility 
along  the  Trans  Alaskan  Pipeline  corridor.  Their  work  was  a  collection  of  man\  other 
investigators’  data  and  results.  Simplifications  of  the  data  presentations  was  made  for  many  of 
the  parameters.  Donovan  and  Singh  (1976)  explain  that  each  of  these  changes,  which  were 
made  so  that  the  procedures  could  be  understood  and  applied  by  geotechnical  engineers  with 
little  or  no  prior  soil  dynamics  experience,  was  compared  with  empirical  data  from  cases  where 
liquefaction  had  been  observed  before  it  was  implemented. 

The  first  criteria  for  decision  making  was  an  examination  of  the  ground  slope  along  the 
route.  Areas  of  slope  less  than  2%  were  not  considered  further  for  liquefaction  potential.  The 
reason  for  this  is  due  to  the  nature  of  the  construction  project.  Bearing  capacity  failures  due  to 
liquefaction  was  not  considered  a  potential  problem.  Bearing  capacity  failures  were  also  not 
considered  problems  of  environmental  concern,  unlike  the  potential  for  flow-slides.  Slopes  in 
excess  of  2%  were  evaluated  with  the  following  non-liquefiable  criteria: 

1.  Soils  consisting  of  clays  and  some  fine  grained  silts  with  sufficient  cohesive  strength  to 
resist  grain  movement  which  could  result  in  the  development  of  excess  pore  pressures. 
This  condition  is  considered  to  be  satisfied  when  the  plasticity  index  number  is  greater 
than  5. 

2.  Soils  are  very  coarse  (large  gravel,  cobbles,  and  boulders)  with  sufficiently  high 
permeability  and  close  proximity  to  drainage  areas  facilitating  rapid  enough  flow  of 
excess  water  to  prevent  build-up  of  pore  pressures. 

3.  Soils  are  not,  or  cannot  become  saturated. 

4.  Frozen  soils  stay  frozen. 

5.  Rock. 

A  flowchart  is  made  up,  and  if  the  soil  still  remains  potentially  liquefiable  after  the  5 
criteria  mentioned,  an  empirical  chart,  relating  dry  density  to  relative  density  by  soil  tyjre  is 
used  to  find  relative  density. 


At  this  point  tluTf  art'  five  sfismic  zones  or  rej^ions  separated  hy  niaxinuim  magnitude 
earthtpiakes  and  expected  duration  i>f  the  ground  shaking.  Each  zone  has  a  loxver  bound 
relatise  density  where  liiiuefaction  potential  is  considered  "high  If  a  density  reading  is 
determined,  design  engineers  are  brought  in  to  further  determine  the  litjuefaction 
susceptibilits'. 

A  Code  of  Aseismic  Design  of  Hydraulic  Structures  for  this  countrs  states  that  ans  silt>'  soil 
which  contains  less  than  15  to  20'<  clay  particles  (the  cla\  particle  is  defined  as  less  than  .005 
mm  in  diameter)  and  has  a  plasticits’  index  greater  than  3,  is  possible  to  reach  liiiuefaction 
during  a  strong  earthquake.  The  code  further  qualifies  this  by  saying  the  water  content  must 
also  be  higher  than  nine-tenths  of  the  .soil’s  liquid  limit. 

The  newest  empirical  method,  and  still  in  its  infant  stage,  is  the  use  of  the  cone 
penetrometer  for  liquefaction  potential  analysis.  .Vlartin  and  Douglas  (1981)  suggest  the 
profiling  of  a  site  using  the  cone  penetrometer  test,  as  opposed  to  the  standard  penetration  te.st 
(SPT).  provides  more  information  fa.ster,  more  reliabh  and  at  a  lower  cost  than  an.v  other 
method  now  available.  Unlike  the  SPT,  where  the  actual  strength  propertv  of  the  soil  being 
measured  is  uncertain  (Seed.  1979),  they  also  insist  that  the  qualitative  determination  of 
strength  variation  throughout  the  site  profile  is  immediately  obtained  from  the  field  records 
due  to  side  friction  measurements  along  the  cone  as  well  as  end  bearing. 

Due  to  the  length  of  existence  of  the  .standard  penetration  test,  and  the  subsequent  broad 
data  base,  re.search  is  now  aimed  at  the  cone  penetrometer  te.st — standard  penetration  te.st 
conversion.  This  emphasis  is  being  focused  on  the  use  of  the  SPT  and  liquefaction  potential 
relationship,  the  methods  of  the  cone  penetration  analysis,  and  the  reliability  of  the  cone 
penetration  (CPT)— SPT  correlations.  Martin  and  Douglas  (1981)  have  found  that  data  .scatter 
in  the  CPT,  in  relation  to  the  SPT,  stems  primarily  from  the  scatter  in  the  SPT  results  or  from 
the  finite  measurement  intervals  of  the  SPT  as  compared  to  the  essentially  continuous  CPT 
measurements.  Even  with  scatter  in  comparative  relations,  the  CPT  method  is  so  rapid  that  a 
large  statistical  data  base  can  be  developed,  allowing  precise  definitions  of  average  site 


characteristics. 


PROBABIJJSTIC  AND  STATISTICAL  METHODS 


Some  of  the  most  reeent  approaches  to  the  problem  contain  concepts  of  statistical  and 
probability  analysis.  (Christian  and  Swiyer  (1975)  presented  a  statistical  approach  that  involves 
the  apparent  relatise  densities  and  ground  accelerations  at  sites  where  licjuefaction  did  and  did 
not  occur  during  an  earth(|iiake.  The  basic  data  was  used  in  a  statistical  analysis  that 
determined  whether  the  soil  would  fit  into  a  li(piefiable  or  non-litpiefiable  category. 
Confidence  levels  or  prohabilits  le\els  were  also  incliitled  within  the  anal\sis  to  determine  how 
good  their  procedure  was.  It  was  pointed  out  that  their  iirobabilitv  levels  were  not  the 
probability  of  liipiefaction.  but  the  level  of  ( oiif idene*'  in  their  dividing  lines  between 
liipiefiable  and  non-li(iuefiahle  sites. 

Yegian  and  Whitman  (1977)  (ireseiited  a  method  of  analvsis  based  on  a  probabilistic 
model.  They  develoiied  a  [larameter  in  their  analvsis  that  is  basicallv  the  ratio  of  the  induced 
cyclic  shearing  stress  to  the  available  strength  of  the  soil.  The  basic  inputs  into  the  parameter 
are  the  magnitude  of  the  earthtpiake  and  the  hypocentral  distance  from  the  site  to  the  causitive 
source.  The\'  include  this  litpiefaction  parameter  in  their  probability  model.  The  iirobabilitv- 
model  gives  the  probability  that  a  site  will  liquefy  under  anv  earthcjuake  loading.  This  is  a 
function  of  the  probability  that  a  site  will  liquefy  given  a  certain  magnitude  earthcpiake  and 
the  probabilitv  of  the  magnitude  earthquake  occurring. 

Youd  and  Perkins  (1977)  have  also  developed  a  procedure  that  is  based  on  probability 
concepts.  The  technirpie  develops  a  liquefaction  potential  map  that  gives  the  relative  potential 
of  a  site  developing  liquefaction.  The  potential  map  is  a  combination  of  two  base  maps,  a 
susceptibility  map  and  an  opportunity  map.  The  susceptibility  map  outlines  the  soil  deposits 
within  a  study  area  that  are  most  likely  to  liquefy.  The  factors  u.sed  to  classify  the  susceptibilitv- 
of  each  deposit  were  the  soil  type  of  the  deposit,  grain  size  distribution,  and  the  age  of  the 
deposit.  A  general  statement  concerning  water  table  depth  was  also  considered  in  their 
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A  liquefaction  opportunity  map  provided  the  seismicity  of  the  study  area.  The  seismicity 
was  determined  from  the  seismic  history  of  the  study  area.  Using  the  seismic  data  and  a 
magnitude-distance  relationship,  a  contour  map  showing  the  return  periods  of  earthquakes 
large  enough  to  cause  liquefaction  was  developed.  The  development  of  the  return  period 
contours  were  ba.sed  on  concepts  from  probability  analysis. 

The  final  potential  map  was  the  combination  of  the  susceptibility  and  opportunitv  maps. 
This  type  of  analysis  is  a  preliminary  guide  for  a  given  study  area.  It  is  not  intended  to  be  used 
as  a  site-specific  analysis  that  could  be  included  in  design  calculations.  It  can.  however,  help  in 
planning  and  site  location  decisions. 

STIFFNESS  METHOD  UTILIZING  IN-SITU  TESTING  PARAMETERS 

Although  the  mechanism  of  liquefaction  has  been  agreed  on.  Dobry  and  Swiger  (1979) 
have  pre.scnted  a  method  of  analysis  which  may  eliminate  the  prexiously  described  problems  of 
laboratory  soil  preparation  for  tests  involving  dynamic  stress  analysis.  An  alternatise  strain 
approach  has  been  proposed  with  the  main  advantages  Ireing  that  the  shear  modulus  at  small 
strains.  G^g^  ~  shear  stress/shear  strain,  can  be  measured  in-situ  b\  geoplnsical  techniques, 
thus  providing  a  direct  measure  of  the  stiffness  of  the  sand  deposit  (G^g.^  =  eV'|.  where  q  = 
ma.ss  density.  Vj  =  shear  wave  velocity).  The  basis  of  this  approach  is  that,  other  things  being 
equal,  .stiffer  sands  are  more  resistant  to  liquefaction. 

As  is  described  in  more  detail  in  a  later  section  of  this  thesis,  a  threshold  strain,  below 
which  no  soil  densification  or  subsequent  pore  pressure  build-up  can  occur,  is  utilized  to 
predict  a  threshold  acceleration  of  a  soil  layer  or  deposit.  Predicted  design  earthquake 
accelerations  can  be  determined  by  computer  analysis,  or  any  other  means,  and  then  compared 
to  the  calculated  threshold  accelerations  to  determine  the  potential  for  pore  pressure  build-up. 
When  threshold  accelerations  are  greater  than  predicted  design  accelerations,  pore  pressure 
build-up  will  occur  and  the  potential  for  liquefaction  should  be  further  investigated  on  a  site- 
specific  basis. 
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OTHER  METHODS 

Two  other  methods  that  do  not  fit  into  the  other  categories  are  presented  b\  Donosan 
(1971)  and  Ghaboussi  and  Dikmen  (1978).  Donovan  (1971)  presented  a  method  referred  to  as  a 
cumulative  damage  approach.  The  method  makes  use  of  .Miner’s  damage  e(piation  and  sums  up 
the  damage  to  the  soil  structure  by  the  cyclic  loading  of  an  earthfjuake.  This  is  analogous  to 
the  fatigue  failure  in  structures.  A  factor  of  safety  is  determined  to  indicate  litpiefaction  or  no 
liquefaction  at  a  particular  site. 

Ghaboussi  and  Dikmen  (1978)  presented  a  procedure  that  models  the  soil  profile  as  a  two- 
phase  flviid-solid  system.  The  method  was  based  on  the  differential  equation  of  motion  of  a 
lumped  ma.ss  system.  The  solution  of  the  equation  of  motion  included  the  nonlinear  properties 
of  the  .soil  and  two  .separate  types  of  damping.  Pore  pressure  distribution  was  monitored  at 
different  depths  by  the  equations  of  motion  to  determine  when  and  where  liquefaction  would 


occur. 
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Chapter  3 

SOIL  SUSCEPTIBILITY  —  THE  QUALITATIVE  STEP 

GENERAL 

A  liquefaction  susceptibility  map  was  developed  for  utilization  in  determining  liquefaction 
potential.  Youd  (1978)  describes  liquefaction  susceptibility  as  the  relative  ease  with  which  the 
material  under  a  particular  site  can  be  liquefied  during  strong  seismic  shaking  and  permanent 
ground  displacement  can  ensue.  The  parameters  considered  in  the  development  of  the 
liquefaction  susceptibility  map  were  the  geotechnically  related  parameters  of  the  unit. 

GEOLOGIC  SETTING 


l,OCATION 

The  China  Lake  Naval  Weapons  Center  is  located  in  two  areas.  The  .southernmost  area. 
Randsburg  Wash,  is  located  south  and  east  of  Searles  Lake  and  Searlcs  Valles.  Due  to  a 
groundwater  level  in  excess  of  40  feet  below  ground  surface,  this  area  was  nrjt  considered  in 
this  study. 

The  second  area,  unofficially  referred  to  as  main  site,  encompasses  some  1000  square 
miles.  The  boundaries  form  a  rectangle,  approximately  42  miles  in  the  north-south  and  26 
miles  in  the  east-west  directions.  The  northern  half  of  the  naval  reservation  is  located  in  the 
southern  extension  of  the  Coso  Range  while  the  eastern  third  of  the  lower  half  reaches  into  the 
Argus  Range.  While  there  may  be  plans  for  the  construction  of  major  facilities  in  these 
mountainous  regions,  the  excessive  depths  to  groundwater  eliminates  any  liquefaction  potential. 

The  remaining  southwest  quarter  of  the  Weapons  Center  is  located  in  Indian  Wells  Valley 
and  extends  west  to  the  base  of  the  Sierra  Nevada  Mountain  Range.  It  is  in  these  lowlands 
where  a  majority  of  the  base's  present  and  future  development  lies.  The  designated  study  area 
for  this  investigation  includes  only  the  southernmost  portion  of  the  described  rectangular  area. 
This  area  covers  approximately  260  square  miles.  (See  Figure  3.1.) 
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GENERAL  GEOLOGY 

Indian  V\'ells  Valley,  including  the  China  Lake  Naval  Weapons  Center,  now  consists  of 
playas  and  alluvial  fans.  During  Pleistocene  this  area  was  occupied  by  perennial  lakes  and 
streams.  The  coalescing  alluvial  fans  extend  from  the  mouths  of  the  surrounding  mountain 
ranges  and  form  gradually  decreasing  alluvial  plains  several  miles  wide.  The  Sierra  Nevadas. 
being  the  major  range  in  the  area,  offer  the  majority  of  the  material  making  up  the.se  fans, 
thus  the  decreasing  slope  is  to  the  east. 

Zbur  (1963)  utilized  geophysical  methods  to  determine  the  thickness  of  the  alluvium  in  the 
central  part  of  the  valley  to  be  nearly  2000  feet.  These  deposits  consist  mainly  of  moderate  to 
well  sorted  sands  with  some  gravels  in  the  higher  regions  of  the  valley,  but  grade  to  silt  and 
.silty  clay  in  the  lower  elevations  where  the  lakes  were  confined  in  more  recent  lime.  Dutcher 
and  Moyle  (1973)  interpret  the  seismic  data  as  revealing  four  main  layers  of  rocks  and  deposits, 
which,  from  deepest  to  shallowest,  have  average  seismic  velocities  of  15,600-16,000  fps;  9.600 
fps;  7,300-7,400  fps;  and  5,700-6,100  fps.  The  authors  describe  these  layers  as  the  highest 
velocity  representing  the  basement  complex,  probably  of  .similar  composition  to  the  Sierra 
Nevada  range,  the  9,600  fps  layer  corresponding  to  the  continental  deposits,  and  the  alluvium 
consisting  of  two  layers,  the  upper  of  which  has  a  velocity  considerably  le.ss  than  the  lower. 
Dutcher  and  Moyle  (1973)  point  out  that  a  refracting  horizon  between  the  deeper  and 
shallower  parts  of  the  alluvium  was  found  consistently  throughout  the  valley  and  proved  to  be 
a  reliable  marker.  The  authors  theorize  that  the  average  velocity  differences  are  caused  by 
compaction  or  consolidation  of  the  deeper  alluvium. 

Dutcher  and  Moyle  (1973)  de.scribe  the  groundwater  basin  as  being  recharged  entirely  by 
internal  surface  drainage.  The  major  .sources  of  this  drainage  are  the  numerous  small 
intermittent  streams  flowing  from  the  canyons  of  the  Sierra  Nevadas.  The  authors  suggest  that 
these  streams  convey  runoff  to  the  fans  where,  after  crossing  the  Sierra  Nevada  fault  zone,  the 
.seepage  recharges  the  groundwater  body.  Additional  recharge  is  generated  from  the  runoff  in 
the  Argus  Range. 
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FAULT  SYSTEMS 

Indian  Wells  Valley,  specifically  the  China  Lake  Naval  Weapons  Center,  lies  in' a  very 
seismically  active  region.  Fault  systems  lying  within  the  valley,  or  nearby,  are  the  main  sources 
of  this  .seismic  activity.  As  described  in  Chapter  5,  background  seismicity  generated  at  a 
di.stance  from  the  study  area  is  also  considered,  but  to  a  lesser  degree.  This  section  will  give  a 
brief  description  of  each  local  fault  considered  for  this  study. 

The  major  source  of  local  seismicity  is  the  Sierra  Nevada  frontal  fault.  This  fault  system 
visibly  extends  a  majority  of  the  length  of  the  Sierra  Nevada  Range.  This  system  can  be  traced 
by  shutter  ridges,  off-.set  streams,  and  other  typical  frontal  fault  geomorphic  features.  Because 
of  its  proximity  to  the  study  area,  only  26  miles  of  the  feature  was  considered  as  a  local  .source 
while  probabilistic  seismicity  generated  by  the  remainder  of  the  fault  was  considered  in  the 
background  seismicity  evaluation. 

The  next  local  source  used  for  this  study  was  the  Little  Lake  fault.  As  described  by 
St.  Amand  (1958),  von  Huene  (1960).  and  Roquemore  (1981),  the  Little  Lake  fault  splass 
eastward  from  the  Sierra  Nevada  frontal  fault  and  continues  south  acro.ss  Indian  Wells  Valley 
to  the  Garlock  fault.  The  trace  of  the  fault  used  in  this  study  was  38  miles  long  and  extends 
directly  through  the  Naval  Weapons  Center  and  this  study  area. 

The  third  local  source  of  .seismicity  used  for  this  study  was  the  Airport  Lake  fault.  Von 
Huene  (1960)  indicates  this  fault  extends  some  18  miles  from  the  north,  where  its  southernmost 
trace  joins  the  Little  Lake  trace  within  the  Naval  Weapons  Center  boundary.  Both  the  Little 
Lake  and  Airport  Lake  faults  are  considered  to  have  experienced  motion  within  the  last  10,000 
years,  thereby  being  classified  as  "active”.  (Roquemore,  1981) 

The  fourth  localized  source  was  the  Argus  fault.  This  sytem  appears  as  a  series  of  splays 
located  on  the  ea.stern  side  of  the  naval  range  area  at  the  base  of  the  Argus  Mountain  range. 
This  splay  region  was  traced  for  20  miles,  encompas.sing  a  wide  band  which  vanishes  ju.st  north 
of  the  actual  study  area. 
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The  fifth,  and  last  local  source  utilized  in  this  study  was  the  Garlock  fault.  This  major 
east-west  trending  fault  extends  approximately  175  miles  from  the  San  Andreas  eastward  to  the 
Death  V'alley  fault  system.  This  feature  is  considered  to  separate  the  Mojase  Desert  and  the 
Basin  and  Range  provinces  (Allen,  et  al,  19651.  Hill  and  Dibblee  (1953)  consider  the  Garlock  to 
be  merely  an  offset  extension  of  one  of  the  major  Transverse  Range  faults.  The  Garlock  trace  is 
located  5  miles  south  of  the  study  area. 


SOIL  TYPES 


GENERAL 

A  .soil  su.sceptibility  map  shows  the  areas  with  soils  that  are  likely  to  liquefy  given  a 
sufficient  magnitude  of  ground  shaking.  As  di.scussed  in  Chapter  2.  .soil  type  is  a  very  important 
parameter  in  analyzing  liquefaction  potential.  A  map  showing  soil  deposits  in  the  Indian  Wells 
Valley  was  prepared  by  Moyle  (1963)  and  was  utilized  in  this  .studs’  to  give  the  information  on 
soil  type  required  for  analysis. 

MOYLE’S  MAP 

Moyle  (1963)  presented  a  detailed  map  of  the  surface  soils  in  the  Indian  Wells  Valles’ 
showing  different  geologic  formations  from  Miocene  to  Holocene.  This  map  provided  the  major 
source  of  soil  data  for  this  .study.  Moyle  is  noted  as  a  particularly  detailed  mapper  and  during 
the  field  checking  for  this  study  there  were  no  discrepancies  noted  in  his  work.  Plate  2  shows 
the  results  of  the  map  that  wa.s  developed  by  Moyle  in  the  study  area. 

Tertiary  volcanics  occur  as  local  hills  and  smaller  outcrops  to  the  northern  regions  of  the 
.study  area.  These  volcanics  extend  from  the  low  lying  foothills  back  into  the  Coso  and  Argus 
ranges. 

Quaternary  deposits  comprise  nearly  all  of  the  valley  floor  region.  Moyle  identified  many 
locations  of  deposits  that  range  in  age  from  Pleistocene  to  Holocene. 


Because  the  age  of  a  deposit  affects  its  liquefaction  characteristics,  the  detailed  mapping  of 
the  Quaternary  units  was  extremely  helpful  in  developing  the  susceptibility  map.  The  seven 
different  .soil  types  that  Moyle  identified  from  the  Quaternary  Period  in  this  study  area  were 
cla.ssified  with  respect  to  liquefaction  potential  and  used  to  produce  the  susceptibilit>  map.  The 
specific  cla.ssification  for  each  soil  type  is  discussed  later. 

Some  of  the  more  highly  susceptible  formations  are  the  .sand  dunes  and  playa  silts  with 
interbedded  sand  layers.  The  dunes  in  the  study  area  are  actively  drifting  and  (p'.ite 
unconsolidated.  The  younger  alluvium  features  a  well  sorted  sand  which  is  slightly  consolidated 
locally,  but  .still  remains  uncemented. 

GEOTECHNICAL  REPORTS 

There  were  only  a  few  geotechnical  reports  at  various  sites  within  the  vallev  that  provided 
information  on  soil  deposits.  These  reports  were  compiled  in  conjunction  with  previous 
foundation  studies  done  in  the  southern  portion  of  the  study  area.  These  studies  were  fairly 
complete  in  their  information,  however,  they  were  few  in  number  and  their  depth  of  studs  was 
shallow,  between  15-20  feet. 

Well  logs,  on  the  other  hand,  were  plentiful  in  number  and  located  over  a  widely 
distributed  area.  The  problems  encountered  with  these,  however,  were  numerous.  First  of  all, 
the  information  on  soil  type  was  very  sketchy  and  suspect.  Further,  the  dates  of  the  loggings 
were  in  the  vicinity  of  1945-1959,  aiid  the  tendency  seemed  to  be  to  lump  a  profile  into  a 
generalized  description.  Also,  the  depths  of  the  borings  were  rounded  off  and  the  personnel 
involved  avoided  detailed  and  consistent  descriptions. 

The  use  of  the  more  complete  geotechnical  reports  is  one  method  that  could  be  very 
helpful  in  refining  classifications  of  susceptibility  to  liquefaction.  Detailed  boring  logs  in  a 
study  area  would  give  more  assurance  in  using  surface  soil  maps  for  the  location  of  liquefaction 
susceptible  deposits.  The  reports  currently  available  are  useful  as  a  supplement  to  a  detailed 
surface  geologic  map,  but  should  not  be  relied  upon  past  their  limited  spatial  accuracy. 


56 


SUSCEPTIBILITY  CLASSIFICATION  SYSTEM 


GENERAL 

The  development  of  the  liquefaction  susceptibility  map  was  based  on  the  factors  discussed 
in  Chapter  2  that  influence  litiuefaction  potential.  The  qualitative  influence  of  several  of  these 
factors  on  liquefaction  and  the  ground  failure  susceptibility  of  various  sedimentary  deposits  has 
been  estimated  b\-  Youd  and  Hoose  (1977)  on  the  basis  of  data  contained  in  the  historical 
records  (Youd  and  Perkins,  1978).  Because  of  the  general  nature  of  the  susceptibilits'  map  not 
all  of  these  influencing  factors  were  considered.  .Vlost  of  the  geotechnical  factors  require  site- 
specific  characteristics  to  determine  their  effect  on  the  liquefaction  potential.  This  prevented 
their  use  on  a  general,  large  scale  map.  There  were,  however,  three  geotechnical  factors  used 
that  were  of  a  nature  that  the  soil  deposits  could  be  clas,sified  as  to  their  relative  susceptibilits 
to  liquefaction.  These  three  factors  were  age  of  deposits,  soil  type,  and  depth  to  groundwater. 
Each  soil  deposit  was  examined  in  terms  of  these  three  factors  and  classified  as  to  whether  it 
had  a  high,  moderate,  low.  or  no  susceptibility  to  liquefaction. 

SOIL  TYPE 

The  first  parameter  considered  was  soil  type  of  the  deposit.  Fine  to  medium  grained  sands 
are  the  most  susceptible  to  liquefaction.  Moyle  did  not  distinguish  between  grain  sizes  within 
sands  and  silts,  such  as  fine,  medium,  and  coarse  grained  material,  but  referred  to  them  only 
as  sand  or  silt  with  a  general  degree  of  sorting.  Therefore,  different  classifications  were  given 
only  on  the  fact  that  the  deposit  was  either  a  .sand,  gravel,  .silt,  or  clay. 

AGE  OF  DEPOSIT 

After  the  soil  type  was  considered,  the  age  of  each  deposit  was  used  to  analyze  the  soirs 
susceptibility.  As  pointed  out  in  Chapter  2,  relatively  recent  deposits  are  the  most  su.sceptible  to 
liquefaction.  The  lower  bound  of  this  study  was  chosen  at  Quaternary  age  deposits  (1.8  million 
years  old;  Pliestocene) .  Deposits  older  than  this  Period  were  classified  as  nonsusceptible  to 


liquefaction.  Younger  depo.sit.s  were  clas.sified  based  on  when  th  'v  were  deposited  within  the 
Quaternary  Period.  Holocene  deposits  (less  than  11,000  years  old)  received  higher  cla.ssifications 
than  did  Pleistocene  deposits  (between  11,000  and  1.8  million  years  old). 

DEPTH  TO  GROUNDWATER 

The  last  geotechnical  factor  that  was  considered  was  the  depth  to  groundwater.  The 
important  influence  of  the  water  level  was  di.sciissed  in  Chapter  2. 

Depth  to  groundwater  data  was  collected  from  the  United  States  Geological  Sursey 
(Lapinski,  1981)  after  a  series  of  wells  were  measured  for  water  levels.  When  data  points  were 
scarce,  a  map,  also  created  by  USGS  (1976),  showing  the  elevation  of  groundwater  above  mean 
sea  level  was  used  in  conjunction  with  a  topographic  map  of  the  same  scale  to  calculate  the 
depth  of  the  water  below  ground  surface.  Latitudinal  and  longitudinal  coordinates  were 
assigned  to  each  data  point  and  used  as  input  into  a  computerized  griding  routine.  From  the 
resulting  grid  system,  a  computer  graphic  program  drew  contour  lines  of  equal  depths  to 
groundwater.  Boundary  conditions  set  for  this  graphic  program  included  intervals  of  2  feet, 
ranging  from  +2.00  to  -40.00  feet.  This  process  is  also  described  in  Chapter  1. 

The  groundwater  contour  map,  as  presented  in  Plate  1,  enabled  three  depth  ranges  to  be 
outlined  and  used  to  further  adjust  deposit  classification.  If  the  depth  to  the  water  level  was 
between  2  feet  above  and  16  feet  below  the  surface  a  higher  classification  was  given  to  a  soil 
deposit  than  if  the  water  level  was  between  17  and  40  feet  below  the  surface.  The  40-foot  level 
was  used  a  lower  bound.  If  the  water  level  is  too  deep  then  the  effective  initial  confining 
pressures  become  so  large  that  liquefaction  is  prevented.  This  was  also  discussed  in  Chapter  2. 
Youd  and  Perkins  (1978)  suggested  that  a  depth  of  33  feet  might  be  a  lower  bound.  Other 
authors  (Shannon,  1966)  (Seed,  1976)  suggest  50  feet  as  a  lower  bound.  A  value  of  40  feet  was 
determined  to  be  a  conservative  yet  realistic  depth. 

INITIAL  SUSCEPTIBILITY 

A  list  of  the  deposits  from  Moyle’s  (1963)  map  is  given  below  with  a  brief  description  of 
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the  soils  contained  in  the  deposits.  A  discussion  on  the  classification  of  liquefaction 
susceptibility  for  each  soil  deposit  is  also  given  in  this  description.  These  classifications  are 
based  onl\’  on  age  of  deposit  and  soil  type.  Depth  to  groundwater  will  be  discussed  later. 

Qof  —  Older  Fan  Deposits.  These  deposits  are  moderately  cemented  to  moderately 
indurated  boulders,  .sand  and  clay,  cemented  with  calcareous  tufa.  On  the  basis  of  soil  typ-" 
and  cementation,  the  su.sceptibility  is  given  a  moderate  cla.ssification.  This  soil  was  deposited 
during  Pleistocene  and,  on  the  basis  of  age,  is  reduced  to  a  low  classification. 

Qol  —  Old  Lacustrine  Deposits.  These  deposits  are  silt,  sand,  and  fre.sh  water  marl, 
cemented  with  calcareous  tufa  cement.  On  the  basis  of  soil  types  and  cementation,  this  deposit 
is  given  a  moderate  cla.ssification  with  liquefiable  lenses.  The  age  of  this  deposit  is  Pleistocene 
which  will  reduce  the  final  classification  to  moderate-low. 

Qls  —  Old  Lakeshore  Deposits.  These  deposits  are  silt,  sand,  clay  and  lime-cement.  On 
the  basis  of  soil  types  and  cementation,  this  deposit  is  given  a  moderate  classification,  with 
some  potentially  liquefiable  lenses.  The  age  of  this  deposit  is  Plei.stocene  which  will  reduce  the 
final  classification  to  moderate-low. 

Qos  —  Old  Dune  Sand.  These  deposits  arc  dunal  sands  of  ni.jderate  consolidation.  On  the 
basis  of  soil  types  and  consolidation,  this  deposit  is  given  a  moderate  classification.  The  age  of 
these  dunes  is  Pleistocene  which  will  reduce  the  final  classification  to  moderate-low. 

Qv  —  Unnamed  Volcanic  Rocks.  These  deposits  are  flows  of  extrusi\e  basalts,  ve.sicular  to 
dense  and  occurring  in  the  higher  regions  of  the  study  area.  The.se  flows  are  Pleistocene  and 
recent.  Volcanic  flows  are  not  susceptible  to  liquefaction. 

Qya  —  Younger  Alluvium.  These  deposits  are  unconsolidated,  moderately  to  well  sorted 
sand,  clay  and  gravel,  and  carried  only  a  moderate  distance  from  the  stturce  area.  On  the  basis 
of  soil  type  and  consolidation,  this  deposit  is  given  a  high  classification.  These  deposits  are 
Holocene,  which  according  to  the  age  maintain  the  classification  of  this  deposit  as  high. 

Qyf  —  Younger  Fan  Deposits.  T'he.se  deposits  are  unconsolidated,  poorly  to  moderately 
sorted  gravel,  sand,  silt  and  mud  flow  debris  locally  derived.  On  the  basis  of  the  presence  of 
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gravel  and  mud  flow  debris  (i.e..  fine  silt  and  clay),  this  deposit  is  given  a  moderate 
classification.  The  age  of  these  deposits  are  Holocene,  thus  maintaining  the  susceptibilits 
classification  of  moderate. 

Qp  —  Playa  Deposits.  These  deposits  arc  unconsolidated  sand,  silt  and  clay  beneath  small 
lake  beds.  On  the  basis  of  soil  types,  i.e.,  large  quantities  of  silt  and  clay  with  localized  sand 
layers,  this  deposit  is  given  a  high  to  moderate  cla.ssification.  The  age  of  these  deposits  is 
Holocene,  thus  maintaining  the  susceptibility  cla.ssification  of  high-moderate. 

Qds  —  Dune  Sand  Deposits.  These  deposits  are  unconsolidated  sand  dunes  which  are 
actively  drifting.  On  the  basis  of  the  soil  type  and  unconsolidation,  this  deposit  is  given  a  high 
cla.ssixication.  The  age  of  this  deposit  is  Holocene,  thereby  maintaining  a  high  classification. 

Qsp  —  Sand  and  Interdune  Playa  Deposits.  These  deposits  are  small  playas  surrounded  b\ 
sand  dunes  which  are  actively  drifting.  Each  of  the.se  depo.sits  have  been  previousK'  described. 
On  the  basis  of  the  .soil  types  and  frequent  occurrences  of  unconsolidated  sands  and  silts,  these 
deposits  are  given  a  high  cla.ssification.  The  age  of  these  deposits  is  Holocene,  thereby 
maintaining  the  high  cla.ssification. 

FINAL  SUSCEPTIBILITY 

Two  of  the  three  geotechnical  parameters  influencing  the  liquefaction  susceptibility  of  a 
soil  deposit  have  now  been  accounted  for  in  the  previous  description.  The  logic  for  each 
decision,  as  suggested  by  Youd  and  Perkins  (1978).  is  illustrated  in  Appendix  1. 

The  third  geotechnical  parameter,  the  depth  to  groundwater,  was  considered  .separately 
from  the  other  two  factors.  This  parameter  depended  on  the  location  of  the  soil  deposit  in 
relation  to  the  depth  to  groundwater  contours.  The  same  .soil  deposit  may  receive  a  high 
susceptibility  classification  in  the  middle  of  the  study  area  where  the  water  is  near  the  surface, 
and  receive  a  moderate  or  low  cla.ssification  in  another  region  where  the  water  level  is 
substantially  deeper. 

The  groundwater  criteria  .set  previously  (0-16',  17-40',  >40')  was  then  combined  with  the 
cla.ssification  results  from  the  age  of  depo.sition  and  the  soil  type  to  form  a  liquefaction 
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susceptibility  table.  Table  3.1  shows  this  susceptibility  matrix.  This  table  illustrates  how  the 
depth  to  groundwater  affected  the  final  classification  of  a  soil  deposit. 


TABLE  3.1.  Susceptibility  Matrix — Estimated  Saseeptibility  of  Soil  Deposits  to 
Liquefaction,  Based  on  Age  of  Deposits,  Soil  Type,  and  Depth  to  Groundwater. 


Depth  to 
Groundwater 

Tertiary 

Quaternary 

Pre-Pleistocene 

Pleistocene 

Holocene 

(feet) 
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EE 

ES 

Qyf 

Qp 
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+2  - 16 

N 

L 

M 

M 

M 

M 

H 

M 

H  M 

H 

H 

17-40 

N 

N 

M/L 

M/L 

M/L 

N 

M 

L 

M  L 

M 

H  M 

>40 

N 

N 

N 

N 

N 

N 

M 

M 

N 

N 

H  -  HIGH  SUSCEPTIBILITY 

H/M  -  HIGH  TO  MODERATE  SUSCEPTIBILITY 

M  -  MODERATE  SUSCEPTIBILITY 

L  -  LOW  SUSCEPTIBILITY 

N  -  NO  SUSCEPTIBILITY 


LIQUEFACTION  SUSCEPTIBILITY  MAP 

By  using  Moyle’s  (1963)  surface  geologic  map  and  transposing  the  computer  generated 
depth  to  groundwater  contour  map  over  it,  the  soil  susceptibilit>’  map  was  laid  out.  At  this 
point  the  .soil  susceptibilits'  matrix  was  used  to  cla.ssify  areas  of  high,  moderate,  and  low 
probable  liquefaction  potential.  The  liquefaction  susceptibility  map  is  presented  in  Plate  2  and 
shows  the  areas  of  relative  liquefaction  susceptibility. 


61 


Chapter  4 

SOIL  SUSCEPTIBILITY  -  THE  QUANTITATIVE  STEP 

GENERAL 

The  (jualitative  step  for  determination  of  soil  susceptibility  to  liquefaction  (Chapter  3) 
oytlined  areas  of  high,  moderate,  and  low  susceptibility.  In  this  chapter,  the  areas  designated 
as  highly  or  moderately  susceptible  are  investigated  further  in  a  more  quantitative  way  to 
formulate  a  liquefaction  potential  map. 

The  procedure  for  developing  this  more  quantitative  analysis  was,  in  part,  discussed  by 
Dobry  and  Swiger  (1979).  The  procedure,  as  Dobry  refers  to  it.  is  called  the  Stiffness  Method. 
As  previously  mentioned  in  Chapter  2,  the  basis  for  this  analysis  is.  other  things  being  equal, 
stiffer  sand  and  silt  deposits  are  less  likely  to  liquefy.  Dobry  (1980)  indicates  that  this  increased 
stiffness  of  the  deposits  can  arise  from  a  combination  of  the  different  parameters  previously 
discus.sed,  such  as  increased  relative  density,  increased  overconsolidation,  stronger  fabric, 
cementation  due  to  geologic  aging  under  pressure,  or  pre.straining  due  to  prior  seismic  history. 

There  have  been  previous  attempts  at  quantifying  liquefaction  investigations  for  the 
purpose  of  learning  the  size  of  an  earthquake  required  to  increase  pore  pressure.  Any  actual 
testing  for  these  results,  unless  purely  empirical,  required  either  standard  penetration  tests  in 
the  field  or  laboratory  tests  in  an  effort  to  determine  the  cyclic  stress  ratio  (rloy)  needed  to 
cause  liquefaction.  As  di.scu.ssed  earlier  in  Chapter  2,  and  as  discus.sed  in  detail  by  Seed  (1979) 
and  Feck  (1979),  the  determination  of  this  stress  ratio  is  riddled  with  uncertainties. 

Dobry  and  Swiger  (1979)  proposed,  ba.sed  on  laboratory  evidence,  that  shear  strains  rather 
than  exclusively  .shear  .stres.ses  control  settlement  and  pore  pressure  build-up  in  silts  and  sands 


during  cyclic  loading.  The\  note  the  soil  parameters  which  increase  stiffness  are  the  same 
which  have  been  shown  to  increase  the  cyclic  strengths  of  sands  under  dynamic  stresses  in  the 
lahorators .  These  detailed  in-situ  parameters,  which  are  nearly  impossible  to  duplicate  in  the 
laborators',  need  not  necessarily  be  known  for  a  particular  site  as  the  stiffness  (shear  modulus, 
or  C»[jiax)  small  strains  can  be  measured  directly  in  the  field  using  geophysical 

techniques. 

The  primars'  information  retpured  for  generation  of  the  (piantitative  susceptibility  map  is 
the  shear  wave  velocits'  of  the  soil  profile,  the  mass  density,  and  the  depth  to  groundwater, 
The  shear  wave  velocity  and  ma.ss  densits'  were  combined  to  determine  the  shear  modulus, 
while  the  depth  to  groundwater  was  used  to  outline  the  potentially  liquefiable  areas  (as 
di.scus.sed  in  Chapter  3)  and  aid  in  the  determination  of  the  total  vertical  stress  due  to 
overburden. 

It  is  not  the  goal  of  this  chapter  to  yield  an  absolute  threshold  ground  surface  acceleration 
capable  of  use  for  design  at  a  particular  .site.  Certain  limitations  and  assumptions  will  be 
implemented  in  order  to  economically  investigate  a  large  geographic  area,  rather  than  a  site- 
by-site  detailed  analysis.  The  end  result  of  this  chapter  will  be  a  grid  system  on  1/2-mile 
centers  indicating  a  threshold  acceleration  capable  of  increasing  pore  pressures  at  a  depth  of 
30  feet  below  the  ground  surface.  The  rea.sons  for  investigating  liquefaction  potential  at  30  feet 
will  be  discussed  in  Chapter  7. 

CONCEPT  OF  THRESHOLD  STRAIN 


Recently,  an  increasing  amount  of  experimental  evidence  has  accumulated  on  the  behavior 

of  .stres.sed  soil  at  the  individual  particle  level.  What  does  soil  undergo  when  strained? 

Figure  4.1  illustrates  the  behavior  of  two  spheres  when  a  normal  load  (N)  and  a  tangential 

load  (T)  are  applied.  The  application  of  the  tangential  load  results  in  a  strain,  y  =  - — —  - . 

Rl  +  R2 

however,  as  long  as  T  <  f  N  (where  f  =  coefficient  of  friction),  the  particles  do  not  slide  past 


each  other.  Mindlin  and  Deresiewiez  (1953)  illustrate  this  load-displacement  behavior  in 


N 


FIGURE  4.1.  Ela.stic  Spheres 
Under  Normal  and  Tangential 
Loads  (From  Mtndlin  and 
Deresiewicz.  1953). 

I'ignre  ‘,.2.  Under  e\(.'licl>  loaded  T.  <iscillating  between  fixed  T*  and  -  T*  (with  T*  <  f  X). 
the  heliasior  is  nonlinear  h> slerefic.  strain-rate  independent,  and  .stable  (Mindlin  and 
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FIGURE  4.2.  Theoretical 
Hysteresis  Loop  Due  to 
O.scillating  Tangential  Force 
for  Two  Spheres  in  Contact 
(From  Mindlin  and 
Deresiewicz.  19-53). 
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Deresiewiez,  1953).  or,  simply  stated,  the  hysteresis  loop  does  not  change  with  the  number  of 
cycles  (Dobry  and  Swiger,  1979).  This  model  is  considered  representative  of  granular  soils  at 
small  strains.  Dobry  and  Swiger  (1979)  suggest  that  the  value  of  the  threshold  strain  for  a  ma.ss 
of  .soil  is  related  to  the  fact  that,  when  the  impo.sed  .strain  equals  the  threshold  strain  (y  =  yj), 
T  =  f  N  and  sliding  occurs  at  a  significant  number  of  grain  contacts.  In  a  similar  graphical 
illustration  of  this,  when  y  >>(,  a  nonlinear,  un.stable  curve  would  result  and  rotate  the  loop  in 
a  clockwise  direction  with  an  increasing  number  of  cycles.  The  physical  result  of  this  model  is  a 
collapse  of  the  soil  .structure  into  a  denser  configuration  until  the  experienced  strain  is  le.ss  than 
the  threshold  strain,  at  which  time  the  hysteresis  loop  will  again  become  stable.  This  behavior, 
for  y  ^yy.  can  best  be  simulated  by  kinematic  strain  hardening.  This  model  suggests.  (1)  y(  is 
independent  of  grain  size,  (2)  yy  increases  with  Oy.  or  N,  (3)  yy  can  be  partialK  dependent  on 
particle  .sorting,  and  (4)  yy  is  quite  dependent  on  the  coefficient  of  friction  and  on  particle 
shape. 

Dobry  and  Swiger  (1979)  related  typical  tests  performed  by  Ortigosa  (1972)  and  Seed  and 
Silver  (1972)  on  dry  sand.  A  shaking  table  apparatus  subjects  a  1-foot-thick  layer  of  sand  to  a 
few  thousand  cycles  of  base  acceleration  and  the  compaction  of  the  sand  layer  is  measured.  It 
has  been  found  that  for  each  relative  density  (Dr)  of  .sand  there  is  a  “critical  acceleration" 
value  (ap),  with  a^  increasing  as  Dr  increases.  If  a  >  a^.,  the  sand  densifies  and  Dr  increases, 
but  the  densification  stops  when  the  condition  a  =  a^.  is  .satisfied. 

Dobry  and  Swiger  (1979)  suggest  that  this  behavior  can  be  explained  by  the  threshold 
strain  concept.  In  these  tests,  the  maximum  cyclic  strain  occurs  at  the  base  of  the  sand  layer 
and  is  equal  to  y  =  (aay)/(gG),  where  a  =  base  acceleration,  g  =  acceleration  of  gravity. 
CTy.  =  vertical  stress,  and  G  =  shear  modulus  of  the  sand  (to  be  explained  in  more  detail 
later).  The  product  (a/g)CTv  is  constant  during  the  test,  but  G  increases  as  the  sand  densifies.  As 
this  is  the  case,  y  decreases  during  the  test  and  densification  stops  when  y  =  yy.  To  restart 
densification,  a  larger  value  of  a/g  is  needed  .so  that  y  >yy  again.  In  this  explanation,  if  a/g  is 
large  enough,  y  is  always  greater  than  yy  and  the  sand  should  keep  densifying  until  it  reaches 


Dr  =  I00‘7c.  This  occurrence  has  been  repf)rfed  b\-  Pyke.  et  al  (1974'  durian  actual  shakinii 
table  tests. 


CONCEPT  OF  MAXIMUM  SHEAR  MODULUS 


The  seismic  shear  strain,  y.  induced  in  the  sand  layer  is  -y  =  r  C.  sshere  G  is  the  shear 
modulus  of  the  soil.  This  value  of  G,  however,  is  not  constant  in  the  soil  throughout  an 
earthquake,  as  it  varies  with  both  the  level  of  strain  and  increasing  pore  water  pressure.  It  can 
be  assumed,  however,  that  for  7  <  -y^  there  is  no  increase  in  the  pore  pre.s.sure  in  the  soil. 
Dobry,  et  al  (1981)  describes  for  the  limiting  ca.se  of  -y  =  "yy  the  corresponding  peak 
acceleration  is  defined  as  the  threshold  acceleration  (ap)f.  The  threshold  strain  and  the 
threshold  acceleration  are  then  related  as  follows: 


7t  = 


(Vt^g 

Gt 


Oy  rd 


(4.1) 


where  Gj  is  the  secant  shear  modulus  of  the  soil  as,sociated  with  the  strain  7^,  and  rd  <  1  is  the 
soil  flexibility  coefficient  defined  by  Seed  and  Idriss  (1971).  This  coefficient  ranges  from 
rd  =  1  at  the  surface  to  rd  =  .92  at  a  depth  of  30  feet. 

This  secant  shear  modulus  is  more  conveniently  expre,s.sed  as  follows: 


Gf  ~  Gjjiax  (G/Gfj,aj()t  (4-2) 

where  G^ax  shear  modulus  of  the  soil  measured  at  very  small  strains  (i.e.,  7  <  10““^ ‘7), 

and  (G/Gjj,ax)  is  the  modulus  reduction  factor  al  the  threshold  strain,  7^. 

Dobry  (1978)  suggests  that  nondestructive  geophysical  methods  for  field  measurements  of 
the  soil’s  shear  wave  velocity  (V^)  induce  very  small  strains  in  the  soil,  approximateh'  the  same 
order  of  those  used  to  define  G,^ax  laboratory.  Therefore,  and  Gj^ax  related  by 

the  basic  wave  propagation  relationship: 

Gmax  ~ 


where  q  is  the  mass  density  of  the  .soil  (total  unit  weight/acceleration  of  gravity). 


DERIVATION  OF  THE  THRESHOLD  ACCELERATION  EQUATION 


Dobry  (1980)  uses  a  simplified  soil  model  of  a  level  sand  site  to  develop  the  proposed 
equation  (see  Figure  4.3).  Seed  and  Idriss  (1971)  and  Seed  (1979)  used  this  same  model  as  part 
of  their  stress  approach  to  liquefaction.  The  static  state  of  stress  on  a  soil  element  at  depth  (z) 


FIGURE  4.3.  Simplified  Soil  Profile  (From  Seed  and 
Indriss,  1971). 


is  defined  by  the  effective  vertical  pressure,  Uy,  by  the  effective  horizontal  pressure,  crj,  = 
kyay,  and  by  the  total  vertical  pressure,  <Ty  =  (7)(Zyy)  +  (7  -  7w)(z  -  z^),  where  7^  = 
62.4  Ibs/ft^  is  the  unit  of  weight  of  water.  The  peak  horizontal  ground  surface  acceleration, 
ap,  and  the  peak  shear  stress,  Zp,  acting  on  the  horizontal  and  vertical  planes  of  the  soil 
element  are  related  by: 


As  discussed  by  Seed  (1979),  a  representative  cyclic  shear  stress  cau.sed  by  the  earthquake, 
Tg,  can  be  estimated  to  be  about  65%  of  the  peak  value,  Tp,  (Dobry,  1980). 

Therefore, 


Tg  =  .65  _  ffy  rd 
g 

and  by  7g  =  Zg/G,  the  associated  cyclic  shear  strain  is: 
g  ^max(^^^max)')'e 

where  G  and  G^iax  de.scribed  in  the  previous  section. 


(4.5) 


(4.6) 
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A  key  experimental  fact  used  in  the  strain  approach  is  the  existence  in  sands  of  a  threshold 
strain  (>().  Several  authors  have  discussed  this  strain  level  (Dobry,  et  al.  1981)  (Park  and  Seed. 
1980)  (Dobry,  1980)  (Powell,  1979)  (Dobry  and  Swiger,  1979)  (Dobrv  and  Grivas.  1978)  (Seed, 
1977)  which  will  be  discussed  in  the  next  .section. 

By  incorporating  this  threshold  strain,  into  equation  4.6,  equation  4.7  becomes: 


(ap)t  a',  rd 

=  .65  — - 

K  ^max  (G/Gmax)T't 


(4.7) 


Rearranging  terms  to  obtain  a  threshold  acceleration  we  get  equation  4.8. 


(^rjlt 


^t  ^max^^^^max^'Vt 
.65  Oy rd 


Dobry  (1980)  de.scribes  equation  4.8  as  the  basis  of  the  .stiffness  method.  If  the  value  of 
Gmax  corresponding  to  the  soil  layer  at  depth  z  is  measured  in  the  field,  then  equation  4.8  can 
be  used  to  determine  the  threshold  peak  ground  surface  acceleration,  (ap)(,  needed  to  start  the 
development  of  excess  pore  pressure  at  that  depth.  A  condition  of  the  design  earthquake 
acceleration,  (ap)  <  (ap)(,  will  indicate  no  pore  pressure  build-up.  The  opposite,  (ap)  ;>(ap)(. 
will  indicate  the  design  earthquake  will  induce  a  pore  pressure  build-up  at  the  indicated  depth 
and  further  liquefaction  studies  are  necessary. 

Dobry,  Yokel,  et  al  (1980)  have  studied  numerous  sites  that  had  liquefied.  Earthcjuake 
peak  accelerations  and  threshold  peak  accelerations  were  calculated  using  existing  standard 
penetration  blow  counts.  For  each  case  studied,  it  was  determined  that  (ap)  ^  (aplj.  The  onK 
two  sites  for  which  (ap)  =  (ap)t  corresponded  to  a  large  magnitude  (M  =  8.3)  earthtjuake 
which  caused  a  very  long  duration  of  shaking  and  eventually  also  liquefied. 


SUSCEPTIBILITY  CLASSIFICATION  SYSTEM 


GENERAL 

The  theory  and  procedure  utilized  for  obtaining  a  peak  threshold  ground  surface 
acceleration  has  been  previously  discus.sed.  The  task  then  was  to  pick  values  to  fulfill 
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equation  4.8  through  field  test  results  and/or  literature  searches  for  pres’iousK  e.stablished 
values  applicable  to  the  conditions  at  the  study  area.  Engineering  judgement  and  a.ssuniptions 
ba.sed  on  the  limited  data  obtained  from  the  field  te.st  program  were  a  fundamental  step  in  the 
application  of  this  method  over  a  large  area.  The  determination  of  each  variable  in 
equation  4.8  is  discmssed  in  this  section  with  the  judgement  and  assumptions  outlined.  The 
variables  determined  are  the  maximum  shear  modulus,  of  the  layer,  the  threshold 

strain,  the  shear  modulus  reduction  coefficient  at  the  threshold  strain,  (G  Gj^ax^tt-  the 
total  vertical  effective  pre.ssure  at  the  layer,  a^, ,  and  the  soil  flexibility  coefficient,  rd. 

DETERMINATION  OF  Gj^ax 

Gfpax  the  maximum  shear  modulus  of  the  soil  at  small  strains,  and  is  defined  as  the 
relationship  of  shear  stress,  r,  divided  by  shear  strain,  y.  This  modulus  is  a  direct  indication  of 
the  stiffness  of  a  soil  layer,  i.e.,  the  layer’s  reaction  in  terms  of  shear  strain  to  applied  shear 
stress.  Glearh'.  the  higher  the  Gp.,3x  value,  the  stiffer  the  soil  profile.  This  stiffness  coefficient 
can  be  measured  in  the  field  by  the  basic  wave  propagation  relationship,  as  discussed  in 

2 

cquati  4. .3  =  Q^s  >• 

In  this  study,  the  shear  ssave  velocity  was  determined  using  the  surface  refraction  method. 
A  12-channel,  signal  enhancing  engineering  geophone  (ES  1200)  apparatus  was  used  for  this 
determination.  The  source  for  the  energv-  impulse  was  a  16-pound  sledge  hammer  striking  an 
end-lwaring  steel  plated  railroad  tie.  A  pickup  truck  was  driven  onto  the  railroad  tie  to 
enhance  the  contact  with  the  ground  (see  Figure  4.4).  As  each  end  of  the  beam  was  struck,  the 
seismograph  recorded  the  seismogram  on  film.  Right  and  left  end  impacts  generated  polarized 
shear  waves  of  a  sinuate  nature  which  were  180°  out  of  phase  of  each  other.  Shear  wave 
velocities  were  chosen  at  points  where  these  sinuous  traces  crossed  each  other. 

Forty-five  sites  were  tested  for  shear  wave  velocity  in  an  effort  to  obtain  an  average 
velocity  for  each  geologic  unit.  The  more  surface  area  a  geologic  unit  covered,  the  more  shots 
were  taken  in  that  unit  for  the  determination  of  this  velocity. 
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FIGURE  4.4.  Generating  Shear  Waves  for  the  Geophysical  Apparatus. 


To  better  compare  individual  localities  within  this  large  study  area,  a  standard  depth  of 
■30  feet  below  the  surface  was  cho.sen  for  analysis  of  susceptibility.  The  rea.sons  for  this  are 
di.scus.sed  in  detail  in  Chapter  7. 

In  the  determination  of  the  shear  wave  velocity  of  the  layer  at  this  depth,  certain 
assumptions  are  made.  First,  it  is  known  that  surface  refraction  is  capable  of  measuring  to  a 
depth  of  1/3  to  1/4  the  length  of  the  string  of  geophones  in  a  uniform,  horizontal  profile 
(Norris,  personal  communications.  1980).  The  string  utilized  in  this  study  was  120  feet  long  (12 
geophones),  thereby  allowing  for  a  .30  to  40  foot  depth  range.  The  assumption  made  at  thi 
point  was  that  there  was  no  "shadow  zone"  in  the  profile.  Briefly,  a  shadow  zone  is  a  layer 
overlaid  by  another  layer  of  higher  velocity.  In  this  condition,  the  shear  velocity  of  the  slower 
deep  layer  will  not  be  detected  by  the  surface  geophones  as  the  waves  will  remain  in  the 


upper,  faster  laser.  The  aceuraes'  of  this  a.ssumption  was  enhanced  by  reviewing  the  available 
geotechnical  boring  logs  in  a  potential  shot  area  for  the  po.ssible  existence  of  such  a  condition. 
When  this  condition  was  noticed,  the  locality  was  not  u.sed  in  the  testing  program. 


With  this  assumption  implemented,  the  velocities  at  30  feet  could  be  measured,  In  the  case 
of  a  constant  velocity  with  depth,  the  calculated  velocity  was  u.sed  for  the  value  at  the  desired 
depth.  In  the  ca.se  of  two  or  more  velocities  measured  in  a  profile,  the  deepest  velocity  (i.e.. 
fastest)  was  assumed  to  be  that  of  the  layer  at  30  feet.  This  assumption  was  verified  by 
calculating  the  depth  of  the  velocity  increase  and  ensuring  that  this  increase  was  not  below  the 
desired  .30-foot  layer.  Equation  4.9  was  used  to  calculate  the  depth  of  the  \  elocity  change  (H). 


D  /V2  -  Vi 

"2  \JV2  + 


(4,9) 


where  D  is  the  “critical  di.stance.  '  or  distance  from  energ\-  source  where  velocities  change 
(from  seismogram),  V2  is  the  faster  velocity,  and  V'|  is  the  slower  velocits'. 

It  is  realized  that  in  reviewing  the  literature  most  authors  agree  that  the  mo.st  popular 
methods  presently  being  u.sed  for  determining  Vj  are,  in  order  of  preference.  (1)  crosshold 
method,  (2)  downhole  method,  and  (3)  surface  refraction  method  (Dobry,  1981)  (Campbell,  et 
al,  1979)  (Stokoe  and  Hoar,  1978)  (Hoar  and  Stokoe,  1977).  However,  due  to  the  large  area 
covered  in  this  study  and  the  potential  astronomical  costs  of  drilling  a  representative  number  of 
boreholes  across  the  Weapons  Center,  the  two  more  popular  methods  will  be  restricted  to 
future  site-specific  investigations.  Table  4.1  illustrates  the  measured  shear  wave  velocities  at  30 
feet  for  the  study  area. 

The  second  component  of  the  equation  4.3  is  mass  deasity  (q).  As  already  di.scussed, 

mass  density  is  defined  as  the  total  unit  weight  divided  by  the  acceleration  of  gravity.  Unit 
weights  were  determined  in  two  ways.  The  most  reliable  method  was,  of  course,  referring  to 
the  existing  geotechnical  logs.  Studies  done  for  foundation  inve.stigations  had  accurate 
measurements  at  depths.  The  second  method  involved  the  u.se  of  the  sand  cone  apparatus  and 
the  existing  well  log  data  (see  Figure  4.5).  The  well  logs,  sketchy  as  they  were,  did  give  an 
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TABLE  4.1.  Measured  Shear  Wave  Velocities  at  a  Depth  of  30  Fec*t. 
Bv  surface  refraction  method. 


- 

_ 

Depo.sit 

V,,  (fps) 

\\  (fpsi 

Deposit 

(fp.s) 

Qol 

1000 

10.37 

Qsp 

1111.11 

888.9 

1162.8 

1222.2 

1428.6 

7,57.6 

Qls 

1100 

IKX) 

1400 

921 

Qos 

1100 

121.3.9 

12.50 

1041.7 

1052.6 

1500 

941.2 

800 

Qya 

1375 

1408.5 

1348.7 

1142.8 

1666.7 

Qds 

1481 

1250 

852 

1315.8 

1333 

808 

1123.6 

1168.8 

947.4 

1463.4 

nil 

1031 

1034.5 

2000 

Qp 

873.8 

915.4 

1081 

880 

833 

909 

indication  of  the  uniformitj'  of  the  profile.  VV'here  it  could  be  determined  that  the  la>er  at  30 
feet  was  approximately  the  same  soil  type  as  that  on  the  surface,  the  sand  cone  method  was 
utilized  for  determination  of  unit  weight.  The  very  rare  condition  involving  a  drastically 
nonuniform  profile  located  in  an  area  without  available  geotechnical  foundation  studies  for 
reference  required  engineering  judgement. 

DETERMINATION  OF 

7f  is  the  threshold  cyclic  shear  strain  discussed  by  many  authors,  as  previously  noted.  Most 
of  these  authors  agree  that  the  key  to  the  .strain  approach  to  liquefaction  analysis  is  the 
development  of  an  accurate  value  for  74.  Richard  Ladd  of  Woodward-Clyde  Consultants,  and 
Ricardo  Dobry  have  worked  on  identifying  this  threshold  strain  (Powell,  1979).  They 
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condiictt'ci  cyclic  friaxial  strain  controlled  tests  on  Monterey  No,  0  sand.  In  these  tests,  the\ 
Siradnally  increased  the  cyclic  strain  of  the  test  until  pore  iiressure  started  to  develop.  Their 
data  clearly  indicated  the  existence  of  a  threshold  strain  value  of  between  1  x  10“^'';  and 
.3  X  lt)“2i;  indeiiendcnt  of  relatise  densitx  (Powell.  1979).  Stoll  and  Kald  (1976)  susmested 
the  value  for  the  threshold  strain  for  saturated  silts  is  sliKhth  lower  (>(  =  ,5  x  ),  In 

rexiewing  the  available  detailed  ijeotechnical  loits,  the  sandy  soil  in  the  studied  layer  contains  a 
moderate  amount  of  silt,  Dtirint;  jiersonal  communications  with  Dobry  (1981),  a  ■>(  = 

.9  X  )()~“‘7  was  sunnested  for  the  describi-d  conditions  in  the  studs  area  ami  he  ajrreed  the 


FIGURE  4  5.  Sand  Cone  Apparatus  for  Mass  Deasity  Determination. 


value  would  he  representative. 


DETERMINATION  OF 

(CJ  has  been  presioush  described  as  the  modulus  reduction  factor  at  the 

threshold  strain  (•>(  =  ),  vshere  G  is  the  st'cant  shear  modulus  of  the  soil  associated  with 

the  threshold  strain  and  is  the  modulus  of  the  soil  measured  at  very  small  strains 

(-y  Recent  laboratory  testing  of  shear  stress  versus  shear  strain  relationships  for 

determination  of  G  values  for  different  strains  has  resulted  in  a  range  of  sariations  for 
G/Gniax-  determined  that  for  =  1  x  (G/G,y^jjx)y(  is  between  .65  and  .95. 

while  for  -yj  =  3  x  ,  (G/G„,3x)7t  between  .45  and  .65  (Dobrv.  et  al.  (981) 

(Iwasaki.  1978)  (Seed  and  Idriss.  1970).  For  the  previously  determined  yy  =  .9  x  for 

this  study,  a  value  of  (C/Cf^^J()7t  ~  selected. 

DETERMINATION  OF  cr; 

Oy  is  defined  as  the  total  effective  vertical  pressure,  or  (y)(z^)  +  (y  -  y^)(z  -  z^.)  = 
Oy ,  where  is  the  depth  to  groundwater,  z  is  the  depth  to  the  soil  element  in  question  (30 
feet),  y  is  the  unit  weight  of  the  soil  above  the  water  table,  y^^.  is  the  unit  weight  of  water,  or 
62.4  lbs/ft3. 

The  unit  weight  of  the  soil  was  determined  as  discussed  in  the  section  regarding 
determination  of  Gjj^^j;.  Depths  to  groundwater  levels  were  read  off  of  the  water  contour  map 
generated  as  discussed  in  Chapter  3. 

DETERMINATION  OF  rd 

As  previous!)'  mentioned,  rd  is  the  soil  flexibility  coefficient  defined  by  Seed  and  Idriss 
(1971).  At  the  ground  surface,  rd  =  1  and  at  a  depth  z  =  30  feet,  rd  =  .92.  It  has  been 
explained  that  shear  stresses  at  any  level  are  due  to  the  upward  propagation  of  shear  waves.  If 
the  soil  column  above  a  given  depth  can  be  considered  to  act  as  a  rigid  body,  the  at 

bedrock  would  correspond  to  the  surface  acceleration  {a)^: 
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’’max  -  Q  ^  (3)s 
rigid 

(force  =  mass  x  acceleration) 

But  the  soil  column  acts  as  a  deformable  body,  therefore 

’’max  “  ’’max 

deformable  rigid 

rd  is  therefore  <  1 

Figure  4,6  contains  the  values  of  rd  Seed  and  Idriss  (1971)  have  determined  for 
of  earthquakes  and  soil  conditions. 
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FIGURE  4.6.  Range  of  Values  of  for  Different  Soil 
Profiles  (From  Seed  and  Idriss,  1971). 
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wide  variety 
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MAP  OF  THRESHOLD  ACCELERATIONS 

By  utilizing  equation  4.8,  with  the  variables  as  described  above,  the  first  state  of  the 
licjuefaction  potential  map  was  prepared  for  the  study  area  showing  the  calculated  peak 
threshold  ground  surface  accelerations  required  to  initiate  excess  pore  pressure  at  a  depth  of  30 
feet.  Again,  due  to  the  limitations  in  available  geophysical  techniques  and  other  assumptions  as 
outlined  previou.sly,  the  acceleration  figures  are  not  intended  to  be  used  in  design  calculations. 
Their  function,  when  combined  with  the  design  ground  accelerations  to  be  discussed  in 
Chapter  5,  is  to  relate  a  relative  susceptibility  as  to  localities  and  approximate  return  intervals 
for  the  conditions  at  the  Naval  Weapons  Center.  The  threshold  acceleration  map  is  presented 
in  Plate  3. 


Chapter  5 

DESIGN  PEAK  EARTHQUAKE  ACCELERATIONS 

GENERAL 

The  second  phase  in  the  deselopment  of  the  liquefaction  potential  map  was  the  inclusion 
of  the  seismicity  of  the  study  area.  Seismic  data  was  included  by  wa\  of  a  computer  program 
called  EQRISK  (McGuire,  1976).  The  purpose  of  this  chapter  is  not  to  write  a  user's  guide  to 
the  program,  or  even  to  attempt  to  explain  the  theory  behind  the -calculations  in  determining 
predicted  accelerations.  The  reader  should  be  able  to  surmi.se  (1)  basicalh’.  what  the  program 
does,  (2)  the  purpose  for  using  EQRISK  in  this  study,  (3)  what  the  input  requirements  for  the 
program  are,  and  (4)  most  importantly,  what  the  output  means  and  how  it  was  used  in  this 
study. 

McGuire  (1976)  indicates  that  the  theory  on  which  this  seismic-risk-analysis  program  is 
based  has  „jen  developed  over  .several  years  (Cornell.  1968,  1971)  (Merz  and  Cornell,  1973).  In 
short,  the  program  takes,  as  input,  a  .series  of  peak  ground  surface  accelerations  to  be  analyzed 
at  each  given  site,  or  grid  point,  as  in  this  study.  The  algorithm  calculates  the  probability  or 
risk  that  the  given  acceleration  will  be  exceeded  by  a  larger  acceleration  at  each  particular  grid 
point.  Once  the  risk  is  calculated,  the  number  of  earthquakes  per  year  for  each  given  fault  that 
could  produce  accelerations  in  excess  of  those  inputed  for  analysis  are  calculated.  The  total 
number  of  expected  earthquakes  at  each  site  producing  the.se  excessive  accelerations  is  the  sum 
of  the  number  of  events  from  all  the  fault  sources  analyzed. 

The  advantages  of  using  EQRISK  are  threefold.  First,  EQRISK  employs  an  easy,  general, 
and  convenient  method  of  specifying  source-area  (fault)  geometry,  which  is  totally  independent 


of  the  site  at  which  seismic  risk  is  to  be  calculated.  The  u.se  of  this  proKram.  in  conjunction 
with  this  study,  negates  the  task  of  having  to  input  the  proximity  of  over  600  grid  points  (sites) 
to  the  five  major  source  areas  (faults)  to  be  discussed  later.  Secondly.  EQRISK  performs  risk 
analysis  at  indis  idual  and  specific  sites,  rather  than  analyzing  regions  or  large  localities.  These 
individual  sites  of  risk  analysis  can  then  be  compared  to  the  already  calculated  threshold 
accelerations  at  the  same  particular  site  in  the  form  of  a  ratio.  {ap)/(ap){  (design  acceleration 
divided  by  thre.shold  acceleration).  The  number  and  location  of  grid  sites  is  totally  up  to  the 
program  user.  Lastly,  EQRISK  calculates  the  actual  number  of  earthquakes  per  year  which 
may  exceed  any  design  accelerations.  This  was  useful  in  this  study  as  it  was  then  very  simple  to 
calculate  recurrence  intervals,  in  years,  of  the  potential  for  liquefaction  at  each  site. 

The  last  two  important  factors  to  be  understood  about  EQRISK  are  the  input  and  the 
output  generated  by  the  program.  The  following  two  sections  will  explain  these  points  in  more 
detail. 


INPUT  REQUIRED 

The  complex  input  requirements  for  EQRISK  are  outlined  in  great  detail  by  McGuire 
(1976),  therefore,  this  .section  will  merely  mention  the  basic  parameters  involved  and  then 
concentrate  on  the  inputed  source  areas,  or  faults.  The  importance  and  complexity  of  these 
remaining  parameters  should  not  be  underestimated,  however,  understanding  the  complete 
program  algorithm  and  functions  of  each  input  variable  were  not  the  intent  of  this  thesis.  The 
input  listing  is  located  in  Appendix  II. 

Another  advantage  of  EQRISK,  though  not  mentioned  previou.sly  because  it  does  not 
directly  effect  the  use  of  the  prograi.i  results,  was  the  ability  to  make  a  simple  and  intelligent 
trade  off  between  efficiency  and  accuracy.  The  variable  NSTEP  was  used  to  establish  an 
“iteration"  nvimber  for  calculating  the  seismic  potential  of  a  source  area  by  its  gross  area. 
Simply  put,  the  higher  NSTEP  will  give  a  more  accurate  spatial  distribution,  and 
correspondingly  a  more  accurate  seismicity,  while  increasing  the  cost  of  the  computer  run  time. 
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The  second  input  parameter  was  a  list  of  “intensities"  to  he  examined.  EQRISK,  being 
strictly  a  probabilistic  algorithm,  has  the  capability  of  examining  Modified  Mercalli  intensits’. 
peak  ground  acceleration,  spectral  velocity,  or  many  other  output  parameters.  The  term 
intensity  is  a  general  term  referring  to  the  form  of  the  desired  output.  The  form  of  the 
intensities  for  this  study,  as  seen  in  the  input  li.sting  in  the  appendix,  was  the  natural  logarithm 
of  the  peak  ground  surface  acceleration.  The  range  of  this  list  encompasses  the  range  of 
threshold  accelerations,  as  calculated  in  Chapter  4.  EQRISK  calculates  the  probabilits'  of 
occurrence  of  some  acceleration  in  excess  of  each  of  these  input  accelerations,  as  will  be 
discussed  further  in  the  next  .section. 

The  input  li.st,  titled  "Risks  Desired,”  .simply  works  backwards  from  the  input 
accelerations  to  be  examined.  This  methodology  assumes  risks  to  be  a  design  parameter.  This 
output  was  not  used  in  this  study. 

The  "Attenuation  Data"  variables  have  the  ability  to  manually  weight  attenuation  due  to 
di.stance  from  the  .source  area  in  the  probabilfstic  calculations.  This  spatial  consideration  was 
automatically  utilized  in  the  calculations,  however,  these  variables  give  the  user  additional 
capabilities  in  simulating  his  specific  study  area  conditions. 

The  input  variables  labled  "No.  of  Gross  Sources"  and  "No.  of  Subsources  in  Gro.ss  Source" 
list  the  number  of  energy  sources  in  the  region  the  user  wants  to  consider  for  potential  peak 
surface  accelerations.  Each  gross  .source,  five  in  this  study,  may  be  divided  into  a  number  of 
rectangular  or  triangular  subsources,  depending  on  the  geometry  of  the  gross  source 
boundaries.  Subsources  aid  in  the  creation  of  a  more  definitive  outline  for  the  gross  source 
area. 

The  next  set  of  input  statements  was  the  most  influential  to  the  programs  output.  As 
mentioned,  five  gro.ss  sources  had  been  outlined  for  this  study.  The  exact  latitudinal  and 
longitudinal  coordinates  of  a  polygon  enclo.sing  each  fault  were  read  into  the  computer.  From 
this  point,  de.scriptive  variables  were  inputed  dictating  specific  characteristics  about  the  fault 
and  how  each  source  was  to  be  treated  .statistically. 
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The  variable  COEF  is  a  coefficient  allowing  subjectise  input  by  the  user  to  modify  a 
source  area  with  known  parameters.  For  instance,  a  coefficient  value  of  -  1  may  be  u.sed  to 
subtract  areas  of  known  low  seismicity  that  happen  to  be  located  in  a  large  quadralateral 
enclosing  a  fault.  A  coefficient  of  1,  such  as  u.sed  in  this  .study,  is  the  variable  waiving  this 
capability. 

The  variable  AMO  and  AMI  are  the  variables  indicating  the  lower-bound  and  upper- 
bound  Richter  magnitude  for  each  gross  source.  These  will  be  identified  later  in  this  section  as 
each  fault  is  discussed. 

The  variable  BETA(d)  is  likened  to  a  characterLstic  "fingerprint"  of  the  seismicity  of  an 
area  in  that  it  describes  the  relative  distribution  of  small  and  large  magnitude  events.  BETA  is 
defined  as  the  natural  logarithm  of  10,  multiplied  by  the  Richter  b  value  for  the  source.  The 
larger  values  of  b  imply  relatively  fewer  large  shocks,  and  visa  versa.  Values  for  b  are 
relatively  constant  for  different  areas,  ranging  from  .67  to  1.29  for  the  continental  United 
States  (Evernda,  1970).  A  value  of  .88  is  cjuite  consistent  with  the  data  of  southern  California 
(Allen  and  others,  1965)  (California  lastitute  of  Technologs'  Seismological  Laboratory  data, 
1981). 

The  variable  RATE  is  the  rate  of  occurrences  of  events  having  magnitudes  greater  than  the 
previou-sly  e.stablished  lower-bound  magnitude  (AMO)  for  each  gross  .source  area.  This  variable 
effectively  describes  the  seismic  activity  rate  of  the  fault.  For  each  gross  source  area,  RATE  is 
in  units  of  the  number  of  events  per  year. 

Parameter  values  to  indicate  “background  seismicity"  may  also  be  input  to  the  program. 
McGuire  (1976)  describes  this  as  seismicity  that  cannot  be  a.ssociated  with  a  specifir^ource 
area.  This  activity-rate  input  is  that  of  a  10,000  km^  area,  i.e.,  an  area  encompassing  a  circle 
around  the  study  area  with  a  radius  of  150  km. 

This  activity-rate,  unlike  the  activity-rates  for  gross  source  areas  being  the  number  of 
events  per  year,  is  in  units  of  number  of  events  per  year  per  10,000  km^.  This  calculation  was 
made  once,  and  the  contribution  to  the  expected  number  of  occurrences  of  peak  ground  surface 
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accelerations  greater  than  those  in  the  list  to  be  examined  was  added  to  the  expected  number 
for  each  site  examined. 


In  order  to  calculate  rate  for  both  individual  faults  and  background  seismicity,  a 
minimum  magnitude  earthquake  to  be  considered  and  a  design  earthquake  mirst  be 
determined.  A  design  earthquake  is  defined  as  the  earthquake  of  highest  credible  magnitude  a 
given  fault  can  produce  ba.sed  on  either  statistics  of  the  length  or  the  displacement  of  faults 
that  have  historical  earthquakes,  or  both  (Roquemore.  1981).  This  data  is  represented  in  the 
form  of  curves  that  plot  the  parameters  of  magnitude  versus  length,  magnitude  versus 
displacement,  as  well  as  other  modifications  (Toucher,  1958)  (lida,  1959,  1965)  (Albee  and 
Smith,  1966)  (Bonilla,  1967,  1970)  (Bonilla  and  Buchanan,  1970). 

Pre.sented  is  a  list  of  the  five  gross  sources,  or  faults  considered  for  this  study,  as  well  as 
the  background  seismicity.  The  data  required  for  the  calculation  of  RATE  is  also  listed. 


Minimum 

Design 

Recurrence 

Gro.ss  .source 

earthquake 

earthquake 

interval 

Reference 

magnitude 

magnitude 

(yrs) 

1, 

Sierra  Nevada . 

2.5 

8.00 

300 

Rociuemore.  personal  comm 

2. 

Little  Lake . 

2.5 

6.80 

1680 

Roquemore,  1981 

3. 

Airport  Lake . 

2.5 

6.56 

2080 

Roquemore.  1981 

4, 

Argus  . 

2.5 

6.50 

1860 

Roquemore.  personal  comm 

5. 

Garlock  . 

2.5 

8.00 

1200 

Roquemore.  personal  comm 

Background  seismicity 

2.5 

8.00 

4138 

Cal.  Tech.,  1981 

From  the  above  data,  the  calculations  for  RATE  start  with  a  relationship  first  presented 
by  Richter  (1958). 

log  njYi  =  a  -  bM  (5.1) 

Simply  put,  ecpiation  5.1  is  an  expression  for  the  number,  niyj.  of  earthquakes  per  year 
having  magnitude  greater  than  M  occurring  in  a  particular  .source  area,  where  a  is  a  constant 
characteristic  of  the  source  area  and  b  was  previously  de.scribed. 

Step  one  was  to  determine  a  value  for  the  constant  a  by  using  what  is  known  about  the 
design  eartht)uake.  For  example.  RATE  was  calculated  for  the  Sierra  Nevada  fault  as  follows; 
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loK  1/300  =  a  -  (.88)(8.0)  (5.1) 

4. .56  =  a 

The  \alue  for  a  will  remain  a  constant  for  the  particular  source  area  while  the  recurrence 
interval  for  the  lower-bound  earthcjuakes  were  evaluated  (RATE): 

log  RATE  =  4.56  -  (.88)(2.5)  (5.1) 

RATE  =  230.61 

Therefore,  it  was  determined  that  the  Sierra  Nevada  fault  would  produce  230.61  earthquakes 
per  year  of  magnitude  greater  than  2.5.  These  calculations  were  repeated  for  each  gro.ss  source 
and  incorporated  into  the  program  as  input. 

OUTPUT  GENERATED 

The  form  of  the  generated  output  is  a  grid  point  by  grid  point  analysis  of  the  probability 
that  an  arbitrary  peak  ground  surface  acceleration  will  exceed  that  acceleration  listed  to  be 
examined  (see  Appendix  III  for  output).  To  exemplify  this,  the  first  grid  point  will  be  discussed 
for  a  listed  acceleration  to  be  examined. 

The  first  site  location  (grid  point),  located  after  the  background  seismicits’  in  the  output  is 

117. .5568°  longitude  and  35.6404°  latitude.  For  illustration,  the  intensity  4.83  will  be 

discu.s.sed.  It  is  this  list  of  intensities  that  is  input  to  be  examined.  The  natural  antilog  of  this 

input  (4.83)  is  calculated  to  be  125.21,  which  represents  a  fraction  of  1  gal  (980  cm/sec^),  or 

125  21 

the  acceleration  of  gravity  Therefore,  this  acceleration  to  be  examined  may  be  more 

commonly  written  as  .13  g.  or  thirteen  hundredths  the  acceleration  of  gravity. 

The  next  line  of  the  output  format  may  be  a  "caution”  statement,  which  is  the  result  of  a 
computer  run-time  cost  versus  efficiency  trade  off  made  by  the  user  in  the  input  mode.  Briefly, 
the  caution  statement  suggests  that  the  value  for  the  variable  NSTEP  does  not  allow  a  sufficient 
number  of  arc  sweeps,  or  "iterations"  to  include  at  least  95%  of  the  total  seismic  gross  area.  If 


the  integration  error  in  the  area  is  greater  than  5%,  the  caution  statement  informs  the  user. 
The  appearance  of  this  statement  merely  reflects  a  slight  reduction  in  the  probabilistic 
accuracy,  but  is  not  a  critical  factor  in  obtaining  reliable  output.  McGuire  (1976)  covers  the 
reason  for  this  statement  and  the  steps  leading  to  it  in  great  detail. 

The  next  line  of  output  is  the  total  number  of  earthquakes  expected  to  be  generated 
annually  which  will  produce  peak  ground  surface  accelerations  equal  to  or  greater  than  the 
listed  accelerations  to  be  examined.  As  discus.sed  earlier.  EQRISK  calculates  risk  or  probability 
associated  with  the  acceleration  to  be  examined.  McGuire  (1976)  explains  that  once  the  risk 
a.ssociated  with  a  given  acceleration  at  a  site  has  been  calculated  for  the  occurrence  of  one 
earthquake  of  arbitrary  magnitude  and  location  in  a  source  area,  the  annual  expected  number 
of  events  from  that  source  area  that  cause  the  given  acceleration  or  greater  is  obtained  b\' 
multiplying  the  .single-event  risk  by  the  expected  number  of  events  during  one  year.  The  total 
expected  number  of  events  causing  the  given  acceleration  or  greater  at  the  site  is  then  obtained 
bv’  summing  the  expected  number  from  each  .source  area.  It  is  this,  the  total  expected  number, 
that  is  reflected  in  the  output. 

The  illustrative  example  used  previously  indicates  .146  E  +  00  (.146)  earthquakes  per  year 
are  expected  with  accelerations  of  at  least  .13  g  (125.21/980).  Another  way  of  stating  this  is 
that  at  least  1  event  of  .13  g  acceleration  will  occur  every  6.85  years  (1/.146  =  6.85). 

The  next  output  line  is  not  used  in  this  study.  This  line  simply  states  the  risk  or  probability 
of  a  random  event  having  a  peak  ground  surface  acceleration  equal  to  or  greater  than  the  one 
.specified  to  be  examined.  In  the  example,  there  is  a  .136  probability  (or  13.6%  chance)  that  an 
event  having  an  acceleration  greater  than  or  equal  to  .13  g  will  occur  each  year. 

The  last  three  lines  are  also  not  used  in  this  study.  The  risks  listed  here  are  “input  risks"  to 
be  examined.  The  intensities  arc  the  logarithms  of  the  peak  ground  surface  accelerations  back 

figured  from  these  given  risks.  Percent  gravity  is  figured  the  same  way  as  was  previously 

(Antilog(Inten.sity)  \ 

- 1  ■ 

980  cm/.sec^  / 
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The  design  peak  earthquake  accelerations  are  added  to  the  threshold  acceleration  map 
(Plate  3).  These  design  accelerations  listed  on  the  map  are  those  of  an  examined  100- year 
interval.  Should  the  reader  desire  an  alternative  interval,  the  same  output  from  EQRISK  ma>' 
be  utilized  as  outlined  previously  in  this  section. 
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Chapter  6 

LIQUEFACTION  POTENTIAL  MAP 

GENERAL 

A  liquefaction  potential  map  is  developed  by  dividing;  the  design  peak  earthquake 
acceleration  for  a  100-year  return  interval,  ap,  by  the  threshold  acceleration  required  for  pore 
pressure  build-up  at  a  depth  of  30  feet,  {ap)(.  The  potential  map,  presented  in  Plate  3,  outlines 
areas  of  relative  potential  for  liquefaction  by  utilizing  this  ratio,  ap/(ap)j,  over  areas  of 
designated  high  groundwater  levels  (less  than  or  equal  to  40  feet  below  ground  surface)  at  the 
Naval  Weapons  Center,  China  Lake,  California. 

Automap  11  (Environmental  Research  Institute,  1981)  was  utilized  to  graphically  illustrate 
the  potential  ratio  by  dividing  the  determined  values  into  11  equal  subdivisions,  as  noted  in  the 
key.  In  general,  the  darker  or  more  detailed  graphics  pattern  corresponds  to  a  higher  potential 
for  liquefaction.  Blank  areas,  as  noted  around  the  perimeter  and  in  two  small  regions  within 
the  potential  zone,  correspond  to  areas  of  low  to  nonexistent  potential.  The  tw'o  blank  regions 
are  labeled  such  due  to  the  frequency  of  clay  and  other  non-liquefiable  soils  as  determined  by 
available  geotechnical  logs.  The  perimeter  region  of  low  potential  is  defined  by  a  low- 
groundwater  level  (greater  than  40  feet  below  ground  surface). 


Chapter  7 
SUMMARY 


DISCUSSION  OF  RESULTS 


GENERAL 

This  study  identified  liquefiable  soil  deposits  at  the  China  Lake  Naval  Weapons  Center 
and  determined  the  relative  potential  these  deposits  have  for  an  increase  in  pore  pressure,  and 
possible  liquefaction  over  a  100-year  design  interval.  This  determination  was  carried  out  in  a 
three  step  approach,  each  step  to  be  discussed  individually  in  this  chapter.  The  locations  of  soil 
deposits  that  arc  su.sceptible  to  liquefaction  are  shown  on  a  licjuefaction  susceptibilits  map. 
These  deposits  are  identified  by  utilizing  chosen  geotechnical  parameters  as  major  influencing 
factors  in  liquefaction  potential.  Deposits  and  areas  expected  to  be  highh  or  moderateb 
susceptible  are  then  subjected  to  a  more  detailed  engineering-type  analysis  for  a  determination 
of  their  threshold  acceleration;  a  minimum  acceleration  recjuired  for  pt)re  pressure  to  begin 
increasing  when  subjected  to  continued  dynamic  loading.  LastK'.  maximum  probab'c 
accelerations,  for  a  100-year  return  interval,  expected  in  these  areas  of  liquefaction 
susceptibility  are  determined  for  comparison  with  the  threshold  accelerations.  The  end  result  is 
a  dimensionless  number,  i.e.,  a  ratio,  indicating  the  anticipated  liquefaction  potential  for  a 
designated  area  for  the  given  return  interval. 

SOIL  SUSCEPTIBILITY  —  THE  QUALITATIVE  STEP 

The  referenced  qualitative  step  yields  a  liquefaction  susceptibility  map  which  is  presented 
in  Plate  2.  The  map  presents  general  areas  where  it  is  highly  probable  that  .soil  conditions  are 
right  for  the  development  of  liquefaction.  A  surface  geologic  map  of  Indian  Wells  Valley, 
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including  Jhe  Naval  Weapons  Center,  prepared  by  Moyle  (1963)  is  irsed  as  a  major  source  of 
reference  in  locating  susceptible  soil  deposits  in  the  study  area.  Alternative  sources  for  geologic 
reference  include  geotechnical  boring  logs  and  well  log  data.  Cross  checks  with  these  sources, 
in  conjunction  with  actual  field  work  performed  by  this  author,  are  u.sed  to  verify  the  accuracs' 
of  the  geologic  information  used  in  this  study.  It  should  be  realized,  however,  that  there  does 
exist  within  these  susceptible  areas  places  or  individual  sites  where  licjuefaction  would  not 
actually  occur.  The  mapping  of  susceptible  regions  on  a  large  scale  did  not  allow  for  the 
prevention  of  misrepresenting  specific  sites  by  indicating  high  classifications  when  the  chance 
for  liquefaction  is  actually  low  or  nonexi.stent.  This  limitation  max  also  appK  to  specific  sites 
located  in  an  area  that  receixed  a  loxv  susceptibility  cla.ssification  xvhen  the  probability  for 
liquefaction  at  the  specific  site  is  quite  high.  This  type  of  misclassification  of  susceptibilitx’  xvas 
unavoidable  and  must  be  realized. 

The  importance  of  the  depth  to  groundxvater  on  liquefaction  potential  xvas  sufficiently 
emphasized.  Because  of  this,  it  must  be  recognized  that  the  accuracy  of  the  computer 
generated  groundxvater  contour  plot  may  be  another  limiting  factor.  Knoxvn  depths  to 
groundxvater  are  verified  on  the  contour  map;  hoxvever,  the  large  areas  of  unknoxvn  data  are 
subject  to  statistical  interpretation  by  the  gridding  and  plotting  routines. 

Upon  completion  of  the  qualitatixe  analysis,  some  general  statements  can  be  made  about 
the  susceptible  areas  at  the  Weapons  Center.  The  vast  majority  of  the  su.sceptible  areas  lie 
within  the  central  portion  of  the  most  recent  China  Lake  bed.  These  areas  of  high  water  table 
are  located  along  G-I  Tower  road  and  Centerline  road.  Although  this  lake-bottom  material 
includes  nonsu.sceptible  soil  such  as  clay,  the  presence  of  silt  and  sand  layers  throughout  the 
depo.sit  creates  a  considerable  risk  for  liquefaction.  Another  region,  which  currently  supports  a 
considerable  amount  of  development,  is  within  a  circle  of  one  mile  radius  around  the  Water 
road  and  Pole  Line  road  intersection.  In  addition  to  playa  deposits,  the  young  alluvium  consists 
of  unconsolidated  sands  and  silts  deposited  in  a  well-sorted  state.  This  material,  in  all 
probability,  has  been  tran.sported  from  the  nearby  Sierra  Nevada  mountains  by  a  variety  of 
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physical  mechanisms,  the  most  predominant  being  water  in  the  form  of  seasonal  runoff  or  flash 
flooding. 

SOIL  SUSCEPTIBILITY  —  THE  QUANTITATIVE  STEP 

The  referenced  quantitative  step  yields  a  map  of  threshold  accelerations  required  to 
initiate  pore  pressure  build-up  at  a  depth  of  30  feet.  The  reason  for  choosing  this  seemingly 
arbitrary  depth  is  twofold.  First,  by  examining  the  relationship  established  in  equation  4.8,  it 
can  be  noted  that  as  the  analysis  increases  in  depth,  the  threshold  acceleration,  (ap)(, 
decreases.  This  may  suggest  deeper  layers  are  more  susceptible  to  liquefaction.  As  previousK’ 
mentioned,  this  is  not  necessarily  true.  Dobry,  et  al  (1981)  also  recognized  this  and  suggests 
that  in  most  sites  the  shear  wave  velocity  also  increases  with  depth.  With  this  in  mind,  the 
plots  of  arrival  times  v'ersus  distances  from  the  source  should  form  a  parabolic  curve,  indicating 
the  velocity  gradually  increases  with  depth.  In  this  study,  this  parabolic  tendency  was  noted  in 
only  2.5-30%  of  the  shots  taken. 

As  .stated  in  Chapter  4,  the  surface  refraction  technique  is  not  widely  utilized  when 
attempting  to  isolate  velocity  characteristics  of  a  profile,  due  to  the  lack  of  control  with  depth. 
It  is  hypothesized  here  that  the  method  employed  is  not  precise  enough  in  its  arrival  times  to 
yield  a  consistent  velocity  increase  with  depth;  rather,  it  yields  an  accurate  approximation  of 
the  average  profile  velocity.  Noted  abrupt  velocity  increases  are  considered  a  change  in  soil 
type.  This  inconsistency  was  also  noted  by  Dobry,  et  al  (1981),  though  to  a  lesser  degree  due  to 
his  use  of  crosshole  and  downhole  geophysical  methods.  Figure  7.1  illustrates  28  in-situ 
measurements  of  shear  wave  velocity  for  sands  and  silty  sands  in  California  by  Campbell,  et  al 
(1979),  Gibbs,  et  al  (1980),  and  by  Dobry.  et  al  (1981).  Twenty  of  the  28  data  jKiints  were 
obtained  using  downhole  techniques;  the  rest  correspond  to  cither  crosshole  or  surface 
refraction  measurements.  The  figure  indicates  .shear  wave  selocitii  s  (or  sand  can  var>  grealK. 
but,  more  importantly,  the  data  scatter  suggests  only  a  relative  increase  in  sehx'it)  -.si(.'  ;'-,.'i. 
and  the  subsequent  potential  for  deviation  from  this  velocity-depth  relationship  It  is. 
therefore,  concluded  that  the  inconsistency,  though  an  annoyance,  is  not  uiuonimon  and 
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SHEAR  WAVE  VELOCITY  IN  FT  S 
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FIGURE  7.1.  Twenty-eight  In-Situ  Measurements  of  Shear  Wave  Velocity  for  Sands  and 
Siltv  Sands  in  California.  (Campbell,  et  al  (1979),  Gibbs,  et  al  (1980),  and  Dobrv,  et  al 
(19k)). 

accuracy  may  be  enhanced  with  a  more  e.xpensive  geophysical  program.  B\  choosing  a  common 
depth  to  all  grid  points  analyzed,  i.e..  30  feet,  overburden  stress  due  to  increasing  depth  at  a 
single  site  would  be  removed  from  the  li.st  of  variables  and  variations  in  threshold  acceleration, 
(ap){.  would  become  a  function  of  the  specific  soil  properties  (mass  density,  shear  wave 
velocity) . 

The  second  reason  for  choosing  this  30-foot  depth  as  a  common  depth  for  analysis  is,  as 
mentioned  in  Chapter  4,  many  authors  suggest  this  as  a  conservative  lower  bound  for  the 
occurrence  of  liquefaction.  Thus,  the  groundwater  table  could  be  at  30  feet  and  surface  failure 
could  .still  be  initiated.  Though  this  study  utilizes  30  feet  in  its  calculations,  it  is  recognized  the 
potential  for  liquefaction  may  extend  to  a  depth  of  40  feet.  This  40-foot  depth  forms  the  actual 
lower  bound  for  the  generation  of  the  liquefaction  potential  map  in  this  study. 

As  mentioned  earlier  in  this  section,  the  threshold  acceleration  values  are  those  required  to 
initiate  pore  pressure  build-up  at  30  feet.  The  map  presents  threshold  acceleration  calculations 
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at  grid  points  on  I/2-mile  centers  throughout  the  study  area.  The  derivation  of  equation  4.8 
has  been  explained  and  the  methods  of  determining  individual  variables  in  the  equation  were 
outlined.  The  shear  modulus,  singled  out  as  being  a  ver>’  important  parameter  and 

was  used  as  the  single  descriptive  variable  for  each  individual  deposit.  The  shear  modulus  was 
defined  as  the  mass  density,  q,  multiplied  by  the  square  of  the  shear  wave  selocitv,  V^. 

The  importance  of  the  depth  to  groundwater  on  liquefaction  potential  is  pointed  out  in 
the  previous  section.  In  addition  to  the  accuracy  limitations  noted,  it  is  important  to  realize 
how  the  presence  of  the  groundwater  affects  the  overburden  stre.ss.  This  effect  is  described  in 
Chapter  2. 

Upon  completion  of  the  threshold  acceleration  calculations,  some  general  statements  can 
be  made  about  the  susceptibility  in  the  study  area.  The  lowest  threshold  accelerations  are 
located  in  areas  of  the  highest  water  table.  Reasons  for  this  are  described  in  Chapter  2  and 
include  the  reduction  in  effective  stress  on  the  examined  30-foot  element.  Soil  deposits  typically 
lower  in  shear  wave  velocity  also  exhibit  lower  threshold  accelerations.  These  "slower"  .soils 
include  the  piaya  (Qp)  and  sand  with  interdune  playa  (Qsp)  deposits.  These  regions  are  in  the 
same  areas  designated  as  particularly  high  risks  by  the  qualitative  method  of  analysis. 

DESIGN  PEAK  EARTHQUAKE  ACCELERATIONS 

The  design  peak  earthquake  accelerations  are  also  illustrated  in  Plate  3.  These  calculations 
represent  the  maximum  expected  acceleration  at  each  grid  point  during  a  100-year  interval. 
These  figures  are  determined  with  the  assistance  of  a  probabilistic  computer  program  called 
EQRISK,  developed  by  Robin  K.  McGuire  (1976).  Individual  variables  used  in  the  program  are 
outlined  in  Chapter  5.  Sources  of  seismic  energy  used  in  the  calculations  are  five  major  faults 
in  the  immediate  area,  as  well  as  a  constant  value  for  background  seismicity  determined  by  the 
California  Institute  of  Technology  Seismic  Laboratory. 

An  obvious  trend  in  the  design  accelerations  is  a  marked  increase  in  expected  magnitude 
progressing  in  a  westward  direction.  The  reason  for  this  trend  is  that  grid  points  closer  to  the 
Sierra  Nevada  frontal  fault  are  exposed  to  a  fault,  not  only  capable  of  a  maximum  magnitude 
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8.00  event,  but  also  quite  active  in  that  it  has  a  substantially  .shorter  recurrence  interval  than 
any  of  the  other  seismic  sources, 

THE  LIQUEFACTION  POTENTIAL  MAP 

The  referenced  liquefaction  potential  map  is  a  graphic  representation  of  the  combining  of 
the  design  peak  ground  surface  earthquake  accelerations,  ap,  for  a  100-year  design  interval 
with  the  peak  threshold  ground  surface  accelerations,  (ap)(,  necessary  for  pore  pressure  build¬ 
up  at  a  depth  of  30  feet.  A  dimensionless  number,  or  ratio,  (ap/(ap)j),  is  assigned  to  each  grid 
point  located  on  1/2-mile  centers  throughout  the  study  area.  Automap  II  (Environmental 
Systems  Research  Inst.,  1981)  is  utilized  as  the  computer  mapping  program  illustrating  zones  of 
increasing  liquefaction  potential  across  the  Naval  Weapons  Center. 

At  this  point,  it  is  very  important  to  clearly  understand  the  physical  meaning  of  the 
“liquefaction  potential  ratio."  This  number  .should  not  be  confused  with  the  standard  safety 
factor  in  that  a  safety  factor  implies  the  occurrence  of  failure  given  a  certain  value,  i.e.,  less 
than  one.  The  suggested  avoidance  of  this  interpretation  is  founded  on  two  factors.  First,  a 
higher  liquefaction  potential  ratio  suggests  a  situation  more  prone  to  ground  failure,  and  thus 
unsafe.  On  the  other  hand,  the  standard  safety  factor,  which  is  utilized  in  many  aspects  of 
engineering,  indicates  failure  potential  in  the  opposite  manner  in  that  a  higher  value  suggests  a 
safer  condition  and  a  value  less  than  one  illustrates  a  driving  force  greater  than  the  resisting 
force  and  ensuing  failure.  The  second  factor,  and  a  very  important  consideration  when 
implementing  the  results  of  this  study,  is  the  realization  that  a  liquefaction  potential  ratio 
greater  than  one  indicates  excess  pore  pressure,  or  a  pore  pressure  build-up  at  a  depth  of  30 
feet  induced  by  the  design  earthquake.  The  possibility  of  liquefaction,  as  it  is  defined  in 
Chapter  3,  and  its  actual  effect  on  engineered  facilities  will  depend  on  factors  such  as  duration 
of  shaking,  the  permeability  and  drainage  boundary  conditions  of  the  liquefiable  layer,  and  the 
dilative  or  contractive  behavior  of  the  soil,  which  in  turn  is  controlled  by  the  relative  density 
and  effective  initial  confining  .stresses  acting  on  the  soil  (Dobry,  et  al,  1981).  Clearly,  the 
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determination  of  damage  potential  on  a  site-specific  basis  will  require  detail  borings  and  a 
layer  b\’  layer  analysis  of  the  effect  the  specific  soil  profile  will  have  on  excess  pore  pressure. 

Upon  completion  of  the  calculations  for  liquefaction  potential  ratios  and  the  liquefaction 
potential  map,  some  general  statements  can  be  niade  about  the  po.ssibility  for  the 
corresponding  hazards  in  the  study  area.  As  also  noted  on  both  the  liquefaction  susceptibility 
map  and  in  comparing  the  threshold  acceleration  values,  the  greatest  potential  for  liquefaction 
occurs  in  the  most  recent  beds  of  the  historic  China  Lake.  Areas  along  Centerline  road  and  G-2 
Tower  road  are  consistently  characterized  by  ratios  greater  than  3.00,  which  indicate  a  strong 
possibility  of  liquefaction.  Other  areas  to  the  west  of  these  roads,  also  located  in  the  more 
recently  active  lake  beds  of  the  valley,  support  ratios  in  excess  of  2.60,  also  a  value  indicative 
of  high  pore  pressures  and  great  potential  for  damage.  Another  area  previously  described  as  a 
po.ssible  hazardous  site,  the  intersection  of  Water  road  and  Pole  Line  road,  show  localized  spots 
of  high  damage  potential  to  the  north  and  east  of  the  intersection,  while  the  regions  to  the  west 
and  south  are  notably  lower  due  to  the  more  competent  soil  type. 

In  light  of  the  creation  of  this  liquefaction  potential  ratio,  the  need  for  numerous  case 
studies  to  help  relate  the  actual  damages  incurred  with  varying  ranges  of  this  number  becomes 
evident.  There  are  simply  not  enough  cases  documented  to  have  an  empirical  feeling  for  the 
damage  a  ratio  of,  for  instance,  2.00  would  actually  inflict  given  a  generalized  soil  profile. 
This  need  has  only  recently  been  recognized  by  Ricardo  Dobry  (Dobry,  et  al,  1981)  and  will 
hopefully  be  recognized  by  other  investigators  in  the  near  future  as  a  valuable  tool  in 
quantifying  liquefaction  potential  in  an  effort  to  reduce  any  future  life  and/or  property  losses. 

RECOMMENDATIONS 

The  development  of  this  liquefaction  potential  map  provides  an  important  insight  into  a 
seismic  hazard  that  exists  at  the  China  Lake  Naval  Weapons  Center.  It  is  recommended  that 
the  potential  map  from  this  project  be  used,  first  of  all,  as  a  general  awareness  study  of  the 
problem  and  an  indicator  that  the  hazard  does  exist,  and  second,  that  there  are  zones  varying 
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greatly  in  their  susceptibility  across  the  study  area.  It  is,  therefore,  considered  that  the  results 
of  this  study  will  be  utilized  in  the  planning  and  development  decisions  made  in  connection 
with  the  growth  of  the  Weapons  Center.  It  is  reiterated  that  any  specific  project  will  require  a 
more  detailed  analysis  for  design,  dependent  on  both  the  magnitude  of  the  structure  and  the 
specific  site  location  with  respect  to  the  predicted  liquefaction  potential  zones. 

A  second  recommendation  is  that  data  for  the  prediction  of  liquefaction  potential  at  the 
Naval  Weapons  Center  be  continually  updated.  Future  geotechnical  boring  logs  collected  for 
any  number  of  reasons  should  be  compiled  and  utilized  for  updating  studies  of  this  type.  Depth 
to  groundwater  data  collected  may  be  used  to  periodically  update  the  computer  contour  map 
generated  for  this  first  study.  Lastly,  as  the  catalog  of  seismic  data  is  expanded  for  the  area,  it 
is  expected  that  the  prediction  of  peak  ground  surface  accelerations  and  their  respective 
recurrence  intervals  will  become  increasingly  more  accurate.  Improvements  in  the 
probabilistic-type  computer  programming  necessary  for  these  predictions  is  also  expected  to 
become  more  reliable  in  the  future. 
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Appendix  I 

SUSCEPTIBILITY  TABLE 


TABLE  A.l.  Estimated  Susceptibility  of  Sedimentaiy  Depmits  to  Liquefaction 
During  Strong  Seismic  Shaking  (From  Voud.  1978). 
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Appendix  II 


DESIGN  GROUND  ACCELERATIONS  FOR  LIQUEFACTION 
POTENTIAL,  CHINA  LAKE  NWC 


NSTEP  =  15  OCALC  =  0  JPBNT  =  0 

List  of  examined  intensities  4,45  4.58  4.70  4.83  4.95  5.08  5,20  5.33  5.45  5. .58  5.70 

Risks  desired  0.2000  0.1500  0.1000  0.0500  0.0100  0.0050  0.00.30  0.0000 

Attenuation  data  =  Cl  C2  C3  SIG.MA  RZERO  RONE  .\AA  BBR 
6.16  0.60  -  1.30  0.00  25.00  0.00  2000.00  0.00 


No.  of  gross  sources  5 

No.  of  subsources  in  gro.ss  sources  3  4  12  7 


Gross  source 

L/S 

COEF 

MO 

Ml 

BETA 

RATE/YR 

EDEPTH 

1 

0 

1.00 

2.50 

8.00 

2.0300 

230.6100 

0.0000 

2 

0 

1.00 

2.. 50 

6.80 

2.0300 

3.6000 

0.0000 

3 

0 

’.00 

2.50 

6.56 

2.0300 

1.8000 

0.0000 

4 

0 

1.00 

2.50 

6.50 

2.0300 

1.7800 

0.0000 

5 

0 

1.00 

2.50 

8.00 

2.0300 

57.6500 

0,0000 

(Background) 

0 

1.00 

3.00 

8.00 

2.0300 

6.0700 

0.0000 

Gro.ss  source 

1 

Subsource  coordinate  data 

117.9418 

35.5825 

117.9108 

35., 5702 

Gross  source 

1 

Subsource  coordinate  data 

117.8846 

35.6786 

117.8487 

35,67.53 

Gross  source 

1 

Subsource  coordinate  data 

117.9269 

35.9298 

117.8920 

35.9371 

Gross  source 

1 

Subsource  coordinate  data 

118.0252 

36.0700 

117.9973 

36.1139 

Gross  source 

2 

Subsource  coordinate  data 

117.4685 

35.5895 

117.4835 

35.6162 

Gross  source 

2 

Subsource  coordinate  data 

117.5639 

35.5610 

117.5623 

35.5898 

Gross  source 

2 

Subsource  coordinate  data 

117.6530 

35.5823 

117.6237 

35.6025 

Gross  source 

2 

Subsource  coordinate  data 

117.7573 

35.7533 

117.7259 

35.7678 

Gross  source 

2 

Subsource  coordinate  data 

117,8633 

35.9273 

117.8317 

35.9402 

Gro.ss  source 

3 

Subsource  coordinate  data 

117.8003 

35.8896 

117.7695 

35.8392 

Gross  source 

3 

Subsource  coordinate  data 

117.7771 

36.1166 

117.7419 

36.1140 

Cross  source 

4 

Subsource  coordinate  data 

117.6856 

35.8748 

117,6151 

35.8740 

Gross  source 

4 

Subsource  coordinate  data 

117.6833 

36.0271 

117.6129 

36.0393 

Gross  source 

4 

Subsource  coordinate  data 

117.7389 

36.1271 

117.6779 

36.1537 

Gross  source 

5 

Subsource  coordinate  data 

119.0255 

34.8238 

118.7753 

34.7735 

Gross  source 

5 

Subsource  coordinate  data 

118.4874 

35.0068 

118.4771 

34.9738 

Gross  source 

5 

Subsource  coordinate  data 

118.0766 

35.2349 

118.0551 

35,2149 

Gross  source 

5 

Subsource  coordinate  data 

117.9854 

35.3590 

117.9513 

35.3100 

Gross  source 

5 

Subsource  coordinate  data 

117.5922 

35.4820 

117.5841 

35.4629 

Gross  source 

5 

Subsource  coordinate  data 

117.0167 

35.6144 

117.0083 

35.5881 

Gross  source 

5 

Subsource  coordinate  data 

116.6572 

35.6413 

116,6653 

35.5753 
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